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Introduction 

General relativity (GR) and quantum mechanics are two of the best verified theories of modern 
physics. While general relativity has been spectacularly successful in explaining the universe 
at astronomical and cosmological scales, quantum mechanics gives an equally coherent physical 
picture on small scales. However, one of the biggest unfulfilled challenges in physics remains to 
incorporate the two theories in the same framework. Ordinary quantum field theories, which 
have managed to describe the three other fundamental forces (electromagnetic, weak and strong), 
have failed for general relativity because it is not perturbatively renormalizable. 

Loop quantum gravity (LQG) [14, 111, 170, 186] is an attempt to construct a mathematically 
rigorous, non-perturbative, background independent formulation of quantum general relativity. 
GR is reformulated in terms of Ashtekar-Barbero variables, namely the densitized triad and 
the Ashtekar connection. The basic classical variables are taken to be the holonomies of the 
connection and the fluxes of the triads and these are then promoted to basic quantum operators. 
The quantization is not the standard Schrodinger quantization but an unitarily inequivalent 
choice known as loop/polymer quantization. The kinematic structure of LQG has been well 



Introduction to Loop Quantum Cosmology 



3 



developed. A robust feature of LQG, not imposed but emergent, is the underlying discreteness 
of space. 

With the aim of obtaining physical implications from LQG, in the last years the application 
of loop quantization techniques to cosmological models has undergone a notable development. 
This field of research is known under the name of loop quantum cosmology (LQC). The models 
analyzed in LQC are mini- and midisuperspace models. These models have Killing vectors which 
reduce the degrees of freedom of full GR. In the case of minisuperspaces, the reduced theories 
have no field-theory degrees of freedom remaining. Although there are field-theory degrees of 
freedom in the midisuperspace models, their number is smaller than in the full theory. Therefore, 
these are simplified systems which provide toy models suitable for studying some aspects of the 
full quantum gravity theory. Moreover, classical solutions are well known (in fact, we are aware 
of very few systems which have closed-form solutions of Einstein equations with no Killing 
vectors) and it is relatively easy to study the effects of the quantization. 

LQC cannot be considered the cosmological sector of LQG because the symmetry reduction 
is carried out before quantizing, and the results so obtained may not be the same if the reduction 
is done after quantization. However by adapting the techniques used in the full theory to the 
symmetry-reduced cosmological models we may hope to capture some of the crucial features of 
the full theory, as well as to obtain hints about how to tackle them. Indeed, one of the generic 
characteristics of LQC is the avoidance of the classical singularity. In the present absence of 
recognized experimental and observational signatures of quantum gravity, this novel and robust 
result has been increasing the hope that LQG may indeed be the correct theory of quantum 
gravity. 

In this article we will review the progress made in the various cosmological models studied in 
LQC in the last few years. A recent review [24] emphasizes aspects that are only briefly men- 
tioned here, such as the "simplified" or "solvable" LQC framework, the details of effective dy- 
namics for FRW models with non-zero curvature and/or cosmological constant, and inflationary 
perturbation theory in LQC. On the other hand, here we focus more on midisuperspaces and 
discuss lattice refinement parametrizations at some length. Before starting, we shall briefly 
recall the main features in the kinematic structures of LQG. Similar ingredients are used in the 
kinematic structure of the LQC models to be discussed later. 

1 Loop quantization 

1.1 Ashtekar— Barbero formalism 

In the Hamiltonian formulation, the four-dimensional spacetime metric is described by a three- 
metric q a b induced in the spatial sections £ that foliate the spacetime manifold, the lapse func- 
tion N and the shift vector N a [1, 160] . Both the lapse N and the shift vector N a are La- 
grange multipliers accompanying the constraints that encoded the general covariance of general 
relativity. These constraints are, respectively, the scalar or Hamiltonian constraint and the 
diffeomorphisms constraint (which is a three- vector) . Therefore, the physically relevant infor- 
mation is encoded in the spatial three- metric and in its canonically conjugate momentum, or 
equivalently, in the extrinsic curvature K a b = C n q a b/2, where n is the unit normal to S and C n 
is the Lie derivative along n [196]. 

LQG is based in a formulation of general relativity as a gauge theory [3, 5, 6, 172, 173], in 
which the phase space is described by a su{2) gauge connection, the Ashtekar-Barbero con- 
nection A l a , and its canonically conjugate momentum, the densitized triad 2 Ef, that plays the 

1 Latin indices from the beginning of the alphabet, a, b, ... , denote spatial indices. 

2 Latin indices from the middle of the alphabet, i,j, . . . are SU(2) indices and label new degrees of freedom 
introduced when passing to the triad formulation. 
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role of an "electric field". To define these objects, first one introduces the co-triad e l a , defined 
as q a b = e l a el5ij, where 5ij stands for the Kronecker delta in three dimensions, and then one 
defines the triad, ef, as its inverse efe J b = SjS^. The densitized triad then reads Ef = y/qef, 
where q stands for the determinant of the spatial three-metric. In turn, the Ashtekar-Barbero 
connection reads [30] A % a = T^ + jK^, where 7 is an arbitrary real and non- vanishing parameter, 
called the Immirzi parameter [123, 124], K l a = K^e^S^ is the extrinsic curvature in triadic 
form, and is the spin connection compatible with the densitized triad. Namely, it verifies 
V i,Ef + eijkT J b E ak = 0, where ey^ is the totally antisymmetric symbol and is the usual spatial 
covariant derivative [196]. The canonical pair (A,E) has the following Poisson bracket: 

{A?(x),Ei(y)} = 87rG 1 6Z6iS(x-y), 

where G is Newton constant and 5{x — y) denotes the three-dimensional Dirac delta distribution 
on the hypersurface E. 

Since the internal Euclidean metric 5ij is invariant under SU(2) rotations, the internal SU(2) 
degrees of freedom are gauge. Therefore, in this formulation of general relativity, besides the 
diffeomorphisms constraint C a and the scalar (or Hamiltonian) constraint C, there is a gauge (or 
Gauss) constraint Qi fixing the rotation freedom that we have just introduced. In the variables 
(Af,E b ), those constraints have the following expression (in vacuum) 3 [186], 

Gi = d a E? + e ijk Y{E ak = 0, 



2\ A Tsm,TSTi\ jpaj jpbk 



Ca - Kb E i - °) 
1 



c = /nrrr^ ^ k ^ " ^ + 7 2 )4n^r] E aj E 
V\ det ( E )\ 

where is the curvature tensor of the Ashtekar-Barbero connection, 
F* ab = d a Ai-d b Ai + e ijk AiA k . 

1.1.1 Holonomy-flux algebra 

The next step is to define the holonomies and fluxes which will later be promoted to basic 
quantum variables. 

The configuration variables chosen are the holonomies of A l a . They are more convenient than 
the connection itself thanks to their properties under gauge transformation. The holonomy of 
the connection A along the edge e is given by 
h e {A) ='p e L dxa Ai(x)T^ 

where V denotes path ordering and Tj are the generators of SU(2), such that [t^Tj] = eij k T k . 

The momentum conjugate to the holonomy is given by the flux of Ef over surfaces S and 
smeared with a su(2)-valued function /*: 



E(SJ) = [ fE?e abc dx b dx c . 
Js 



The description of the phase space in terms of holonomies and fluxes is not only suitable for 
its transformation properties, but also because these objects are diffeomorphism invariant and 
their definition is background independent. Moreover, their Poisson bracket is divergence- free 

{E(S, /), h e (A)} = 27rG 7 e(e, S)f i Tih e (A), 

where e(e, 5) represents the regularization of the Dirac delta: it vanishes if e does not intersect S, 
as well as if e C S, and |e(e, S)\ = 1 if e and S intersect in one point, the sign depending on the 
relative orientation between e and S [14]. 



3 In the presence of matter coupled to the geometry, there is a matter term contributing to each constraint. 
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1.2 Kinematic Hilbert space 

In LQG, the holonomy-flux algebra is represented over a kinematical Hilbert space that is 
different from the more familiar Schrodinger-type Hilbert space. It is given by the completion of 
the space of cylindrical functions (defined on the space of generalized connections) with respect 
to the so-called Ashtekar-Lewandowski measure [13, 15, 17, 27]. We give a very brief description 
of this kinematical Hilbert space below, while the details can be found in [14, 111, 170, 186] 
(and references therein). 

A generalized connection h e (A) = A e is an assignment of A S SU(2) to any analytic path 
e C S. A graph T is a collection of analytic paths e C S meeting at most at their endpoints. 
We will consider only closed graphs. The point at which two edges meet is called a vertex. Let 
n be the number of edges in T. A function cylindrical with respect to V is given by 

MA) :=fr(A ei ,...,A en ), 

where f-p is a smooth function on SU(2) n . The space of states cylindrical with respect to T are 
denoted by Cylr- The space of all functions cylindrical with respect to some r G X is denoted 
by Cyl and is given by 



Cyl = UCyl r . 



Given a cylindrical function ipr(A) G Cyl, the Ashtekar-Lewandowski measure, denoted 
by /io, is defined by 



/ diM>W{A)\ ■= [ Y[dh e fr{A ei ,...,A en ), VVr(A), 



'SC/(2)« ecr 

where dh is the normalized Haar measure on SU (2). Using this measure we can define an inner 
product on Cyl: 



(ipr^'r) ■= (fr[A ei ,...,A en ),g r ,(A ei ,...,A em )) 

= I [J dh£fr ( Ae i ,---,A en )g r ,(A ei ,..., A em ) , 

JSU(2) n - 



isu(2y 

where Hpp' is any graph such that T C Hpr' and V C Hrr'- Then, the kinematical Hilbert space 
of LQG is the Cauchy completion of Cyl in the Ashtekar-Lewandowski norm: %kin = L 2 (A, cfyzo). 

A basis on this Hilbert space is provided by spin network states, which are constructed as 
follows. Given a graph T, each edge e is colored by a non-trivial irreducible representation -Kj e 
of SU (2). Spin network states are cylindrical functions with respect to this colored graph. They 
are denoted by T s := T r j(A) where j = {j e }- Then, every cylindrical function can be expanded 
in the basis of spin network states. 

On Cylp the operators representing the corresponding holonomies act by multiplication, while 
the operator representing the flux is given by 



E r (S,f) = i27rGh^2e(e,S)Tr^f n A e ^j . 



ecr 

To obtain the quantum version of the more general operators, they have to be first rewritten 
in terms of the basic holonomy-flux operators. Note that the quantum configuration space is not 
the space of smooth connections but rather the space of holonomies (or generalized connections). 
Since the Ashtekar-Lewandowski measure is discontinuous in the connection, there is no well- 
defined operator for the connection on %kin- Consequently, the curvature must be defined in 
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terms of holonomies before it can be promoted to a quantum operator. The strategy in the full 
theory is to define any general quantum operator via regularization as follows (see [14, 186] for 
details): 

• the spatial manifold S is triangulated into elementary tetrahedra; 

• the integral over S is replaced by a Riemann sum over the cells; 

• for each cell, we define a regularized expression in terms of the basic operators, such that 
we get the correct classical expression in the limit the cell is shrunk to zero; 

• this is promoted to a quantum operator provided it is densely defined on 'Hkin- 

In the subsequent sections we shall see how the same strategy is applied for defining the quantum 
operators in LQC. One significant difference is that in the full theory the final expressions are 
independent of the regularization, while in the symmetry-reduced models the regularization (i.e., 
the size of the cells) cannot be removed and has to be treated as an ambiguity. However, we 
can fix the form of the ambiguity by taking hints from the full theory. 

One of the most interesting features of LQG is that the spectra of the operators representing 
geometrical quantities like area and volume are discrete. Discrete eigenvalues imply that the 
underlying spatial manifold is also discrete at least when we are close to the quantum gravity 
scale. This is a feature of the quantization scheme and it also plays and important role in the 
singularity avoidance in LQC minisuperspace models. 

This is the kinematical structure of LQG. However we are interested in physical states, i.e. 
states which are annihilated by the all the constraints. To obtain the physical Hilbert space 
we now need to solve the quantum constraints. The Gauss constraint is easy to solve and we 
can obtain gauge invariant Hilbert space spanned by the gauge invariant spin networks. The 
infinitesimal diffeomorphism constraint cannot be expressed as a self-adjoint operator on %kin- 
However we can consider finite diffeomorphisms and the solutions to the finite diffeomorphism 
constraint are obtained via group averaging. It turns out that these solutions do not lie in %kin 
but in Cyl*, the algebraic dual of Cyl. 

In the construction of the Hamiltonian constraint operator we face a number of problems 
(see [72] and references therein for details). Although a well-defined Hamiltonian constraint ope- 
rator can be constructed which satisfies an on-shell anomaly-free quantum constraint algebra, 
the quantization procedure suffers from a number of ambiguities: in the choice of the regulators, 
in the transcription in terms of basic quantum variables, and in the choice of curvature appro- 
ximants. Also the domain of the Hamiltonian constraint operator is not known. Efforts have 
been made to reduce the ambiguities by studying the off-shell closure of the constraint algebra 
and by trying to find the correct semiclassical limit, but no significant progress has been made 
so far. So, although we have a well-defined full quantum theory of gravity at the kinematical 
level, the physical Hilbert-space construction is beset by a number of open problems and is not 
yet complete. 

LQC tries to study some of the features of Loop quantization while avoiding the problems 
of the full theory. As we shall see later, the programme of LQC tries to closely follow the 
same steps, as far as possible, in the much simpler case of cosmological models with no (or at 
most one) field-theory degrees of freedom. In minisuperspace models it is possible to go beyond 
the kinematics and construct the physical Hilbert space. Another useful procedure developed 
to study the effect of the underlying discreteness is the use of effective equations to study 
homogeneous cosmologies and perturbations therein. This has opened up a large number of 
systems to semiclassical analyses. It is hoped that lessons learned from LQC can give hints 
about how to tackle the issues being faced in LQG. 
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2 Plan of the review 

Significant progress has been made in the study of a number of cosmologies in LQC. Here, we 
shall give an overall account of various facets of LQC, outlining technical aspects, reviewing 
the results achieved and indicating the directions of further research. The rest of the paper is 
divided into three parts. 

In Part I we discuss LQC minisuperspace models. The simplest cases of minisuperspace 
are Friedmann-Robertson-Walker (FRW) models, which are homogeneous and isotropic. The 
kinematical quantization programme followed for these models will be discussed in detail, using 
the example of flat FRW. We also describe the results obtained in the physical Hilbert space 
including the dynamical singularity resolution and the bounce. Open and closed FRW models, 
with and without a cosmological constant, are briefly discussed. The next level of complication, 
Bianchi models, consists in removing the assumption of isotropy. In this case, our illustrative 
example will be the Bianchi I model but we also indicate the work done so far for Bianchi II 
and Bianchi IX cases. 

Then, Part II focusses on the LQC of midisuperspace models which are neither homogeneous 
nor isotropic. We describe the only case whose loop quantization has been studied in some 
detail, the linearly polarized Gowdy T 3 model. Two contrasting approaches have been taken in 
the study of this model. In the first approach, the degrees of freedom have been separated into 
homogeneous and inhomogeneous sectors. The homogeneous sector is quantized using the tools 
developed in LQC, while the inhomogeneous sector is Fock quantized. In the second approach, 
the model is studied as a whole mimicking the steps of LQG. We describe and compare both 
procedures. 

Finally, in Part III we discuss the programme of effective dynamics developed in LQC. In 
contrast to the previous two parts, this approach aims to incorporate the effects of the discrete 
geometry as corrections to the classical equations. In this way it may be possible to link LQC 
to phenomenological evidence. 

In the end we summarize the current directions of ongoing research. This review is intended 
as an introduction of the main results achieved in the field in the past few years, especially in 
the Hamiltonian formalism, and it does not cover more recent work being done in the area of 
cosmological perturbations, phenomenology, and spin-foam cosmology. We will comment about 
these and other lines of research in Sections 9 and 10. 

Part I 

Minisuperspaces in loop quantum cosmology 

LQC [2, 4, 44, 149] adapts the techniques developed in loop quantum gravity [14, 170, 186] to the 
quantization of simpler models than the full theory, as minisuperspace models. Minisuperspace 
models are solutions of Einstein's equations with a high degree of symmetry, so much so that 
there are no field theory degrees of freedom remaining. They lead to homogeneous cosmological 
solutions all of which suffer from a singularity where the classical equations of motion break 
down. Since, after quantization, these are essentially quantum mechanical systems, they serve 
as good toy models for testing the predictions of LQG. 

In LQC, we start from the classically symmetry-reduced phase space and then try to apply 
the steps followed in LQG to these systems. Owing to simplifications due to classical symmetry 
reduction, many technical complications typical of LQG can be avoided, and the quantization 
programme can be carried out beyond what has been achieved so far in the full theory. The fact 
that there is a well-defined full theory which tells us that the underlying spatial geometry is 
discrete is a crucial ingredient in the formulation of LQC. A significant achievement of LQC is 
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the development of a well-defined quantum theory for cosmological models where the classical 
singularity is absent. This resolution of the classical singularity is a robust feature of LQC as it 
is seen in all the minisuperspace models studied so far, as well as under various choices made 
in addressing the ambiguities arising in quantization. In this part we shall review the LQC of 
various known minisuperspace cosmological scenarios. 

3 Friedmann— Robertson— Walker models 

LQC started with the pioneering works by Bojowald [38, 46, 47, 48, 49], that showed the first 
attempts of implementing the methods of LQG to the quantization of the simplest cosmological 
model: the flat Friedmann-Robertson-Walker (FRW) model (homogeneous and isotropic with 
flat spatial sections), whose geometry is described by a single degree of freedom, the scale 
factor. This system, even if very simple, is physically interesting since, at large scales, our 
universe is approximately homogeneous and isotropic. In addition, cosmological observations 
are compatible with a spatially flat geometry. 

After the early papers by Bojowald, the kinematic structure of LQC was revised and more 
rigorously established [9], which made it possible to complete the quantization of the model 
in presence of a homogeneous massless scalar field minimally coupled to the geometry, as well 
as to study the resulting quantum evolution [12, 19, 20, 21]. Classically, this model represents 
expanding universes with an initial big bang singularity, where certain physical observables, 
such as the matter density, diverge. Remarkably, the quantum dynamics resolves the singularity 
replacing it with a quantum bounce, while for semiclassical states it agrees with the classical 
dynamics far away form the singularity. Therefore, even though this is the simplest cosmological 
model, its loop quantization, also called polymeric quantization, already leads to relevant results, 
the most important one being the avoidance of the singularity. 

Using the example of the flat FRW model coupled to a massless scalar, we shall discuss 
in detail the basics and the mathematical structure of LQC, adopting the so-called improved 
dynamics prescription [21]. 

3.1 Classical phase space description 
3.1.1 Ashtekar Barbero formalism 

The classical phase space in the presence of homogeneity is much simpler than the general situ- 
ation described in the introduction. In homogeneous cosmology, the gauge and diffeomorphisms 
constraints are trivially satisfied, the Hamiltonian constraint being the only survivor in the 
model. Moreover, for flat FRW the spin connection vanishes. In this case, the geometry part of 
the scalar constraint in its integral version is Cg Tav (N) = iVC grav , with 



Since flat FRW spatial sections £ are non-compact, and the variables that describe it are 
spatially homogeneous, integrals such as (3.1) diverge. To avoid that, one usually restricts 
the analysis to a finite cell V. Owing to homogeneity, the study of this cell reproduces what 
happens in the whole universe. When imposing also isotropy, the connection and the triad can 
be described (in a convenient gauge) by a single parameter c and p, respectively, in the form [9] 




'bk 






4 The lapse function TV goes out of the integral due to the homogeneity. 
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Here we have introduced a fiducial co-triad °e l a that we will choose to be diagonal, °e % a = 5 l a , 
and the determinant y/°q of the corresponding fiducial metric. The results do not depend on 
the fiducial choice. With the above definitions, the symplectic structure is defined via, 

, , 8vrG7 
{c,p\ = — ^— ■ 

The variable p is related to the scale factor a commonly employed in geometrodynamics through 
the expression a(t) = \/\p{t)\V 1 ^ 3 . Note that p is positive (negative) if physical and fiducial 
triads have the same (opposite) orientation. 

On the other hand, a (homogeneous) massless scalar field (ft, together with its momentum P^, 
provide the canonical pair describing the matter content, with Poisson bracket {(ft, P^} = 1. 
Then, the total Hamiltonian constraint contains a matter contribution beside the geometry one, 
given in equation (3.1), and reads 

C = C grav + C mat = --c 2 y / \i\+ 8vrG-^ = 0, (3.2) 
where V = |p| 3//2 is the physical volume of the cell V. 



3.1.2 Holonomy-flux algebra 

When defining holonomies and fluxes in LQC, and in the particular case of isotropic FRW 
models, owing to the homogeneity it is sufficient to consider straight edges oriented along the 

1/3 

fiducial directions, and with oriented length equal to fiV , where \i is an arbitrary real number. 
Therefore, the holonomy along one such edge, in the i-th direction, is given by 

hl{c) = = cos (f) l + 2sm(^)n. 

Then, the gravitational part of the configuration algebra is the algebra generated by the matrix 
elements of the holonomies, namely, the algebra of quasi-periodic functions of c, that are the 
complex exponentials 

A/;( C ) = e^ c . 

In analogy with the terminology employed in LQG [14, 186], the vector space of these quasi- 
periodic functions is called the space of cylindrical functions defined over symmetric connections, 
and it is denoted by Cyl s . 
In turn, the flux is given by 

E{SJ)=pV- 2 ' z A SJ , 

where As f is the fiducial area of S times an orientation factor (that depends on /). Then, the 
flux is essentially described by p. 

In summary, in isotropic and homogeneous LQC the phase space is described by the variab- 
les A/^(c) and p, whose Poisson bracket is 



{Af^c),p} = 
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3.2 Kinematical structure 

Mimicking the quantization implemented in LQG, in LQC we adopt a representation of the alge- 
bra generated by the phase space variables M^{c) and p that is not continuous in the connection, 
and therefore there is no operator representing c [9]. More concretely, the quantum configuration 
space is the Bohr compactification of the real line, MBohn an d the corresponding Haar measure 
that characterizes the kinematical Hilbert space is the so-called Bohr measure [192]. It is simpler 
to work in momentum representation. In fact, such Hilbert space is isomorphic to the space of 
functions of fx G M that are square summable with respect to the discrete measure [192], known 
as polymeric space. In other words, employing the kets |^) to denote the quantum states A/^(c), 
whose linear span is the space Cyl s (dense in RBohr)) the kinematical Hilbert space is the com- 
pletion of Cylg with respect to the inner product = S^i. We will denote this Hilbert 
space by 7i gra , v . Note that "Hgrav is non-separable, since the states \fi) form a non-countable 
orthogonal basis. 

Obviously, the action of on the basis states is 

•A^'l^) = |A* + A*')- 
On the other hand, the Dirac rule [N^p] = ih{Af/j,(c),p} implies that 

pIm) = p(m)Im)> p(p) = 47r/ 3 P ' 7 /^ 

where Zpi = \/Gh is the Planck length. As we see, the spectrum of this operator is discrete, as 
a consequence of the representation not being continuous in ji. Due to this lack of continuity, the 
Stone-von Neumann theorem about the uniqueness of the representation in quantum mecha- 
nics [180, 195] is not applicable in this context. Therefore, the loop quantization of this model is 
inequivalent to the standard Wheeler-DeWitt (WDW) quantization [99, 199], where operators 
have a typical Schrodinger-like representation. In fact, while the WDW quantization fails in 
solving the problem of the big bang singularity, the loop quantization is singularity free [20, 21], 
as we will see later. 

For the matter field, we adopt a standard Schrodinger-like representation, with (j> acting by 
multiplication and Pa, = —ihd^ as derivative, being both operators defined on the Hilbert space 
L 2 (R, d(f>). As domain, we take the Schwartz space <S(M) of rapidly decreasing functions, which 
is dense in L 2 (M, d(j)). The total kinematical Hilbert space is then %kin = %grav <8> L 2 (R, d(j)). 5 

3.3 Hamiltonian constraint operator 

3.3.1 Curvature operator and improved dynamics 

Since the connection is not well defined in the quantum theory, the classical expression of the 
Hamiltonian constraint, given in equation (3.2), cannot be promoted directly to an operator. In 
order to obtain the quantum analogue of the gravitational part, we follow the procedure adopted 
in the full theory. We start from the general expression (3.1) and express the curvature tensor 
in terms of the holonomies, which do have a well-defined quantum counterpart. 
Following LQG, we take a closed square loop with holonomy 

= Kh^r\^)-\ 

5 Note that the basic operators defined above are in the tensor product of both sectors (geometry and matter), 
acting as the identity in the sector where they do not have dependence. For instance, the operator p defined on 
Cyl s Cg> 5(K) really means p (g> 1. Nonetheless, for the sake of simplicity we will ignore the tensor product by the 
identity. 
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that encloses a fiducial area Aq = ji 2 V . The curvature tensor then reads [9] 



Fl b = -2 lim tr -r l °ei°e k b . 



(3.3) 



This limit is classically well defined. However, in the quantum theory we cannot contract the 
area to zero because that limit does not converge. 

Since we have a well defined full theory (unlike WDW quantization), we can appeal to the 
discretization of geometry coming from it. In LQG, geometric area has a discrete spectrum with 
a non- vanishing minimum eigenvalue A [16, 171]. This suggests that we should not take the null 
area limit, but consider only areas larger than A. Then, we contract the area of the loop till 
a minimum value ^4n min = £l 2 Vo^, such that the geometric area corresponding to this fiducial 
area, given by the flux E(n m i n , / = 1) = pp?, is equal to A. In short, the curvature is defined 
by the regularized expression 



This choice of fl is usually called improved dynamics in the LQC literature [21]. Note that the 
smaller the value of fl is, or equivalently the bigger the value of \p\ is, the better equation (3.4) 
approximates the classical expression (3.3), so that both expressions agree in the regime in which 
the area of the cell under study is large enough. Finally, the curvature operator is obtained by 
promoting equation (3.4) to an operator. Let us remark that there are two kinds of ambiguities 
in the definition of this operator. On the one hand, the value of the parameter fl, that is fixed 
by the improved dynamics prescription, as we have just explained. On the other hand, we also 
have the ambiguity in the SU(2) representation we use for calculating the trace. As usual in 
LQC [44], we will compute the holonomies in the fundamental representation of spin 1/2. 

Note that terms of the kind J\fp_ = e % ^ c / 2 contribute to h^,, . In order to define the operator 

A/a = e*^ c / 2 , it is assumed that this operator generates unit translations over the affine param- 
eter associated with the vector field /2[p(/z)]d M [21]. In other words, we introduce a canonical 
transformation in the geometry sector of the phase space, such that it is described by the variable 
b = Hfic/2 and its canonically conjugate variable v(p) = (2-7r7/p 1 v / A) _1 sgn(p)|p| 3//2 (sgn denotes 
the sign), with {b, v} = 1. The variable v(p) = v\p(fi)] indeed verifies d v = fl(fj,)d^. Then, we 
relabel the basis states of % gra v with this new parameter v that, unlike fi, is adapted to the 
action of Jvn. In fact, introducing the operator ii with action v\v) = v\v), it is straightforward 

to show that Afp\v) = \v + 1), so that the Dirac rule [e lb / fi , v]\v) = ih{e tb / h : v}\v) is satisfied. On 
the other hand, we obtain p\v) = (27T7/p 1 \/A) 2 ^ 3 sgn(f )\v \ 2 / 3 \v). 

It is worth mentioning that the parameter v has a geometrical interpretation: its absolute 
value is proportional to the physical volume of the cell V, given by 

—-3/2 * „ ,— 

V =\p\ , V\v) = 2^l\ x ^\v\\v). 

The quantization within the prescription (3.5) meant an important improvement for LQC [21]. 
Earlier, it was assumed that the minimum fiducial length was just some constant \i a related 




(3.4) 



where fj,, characterizing the minimum area of the loop, is given by the Ansatz 




(3.5) 
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to A [9]. However, the resulting quantum dynamics was not successful, inasmuch as the quan- 
tum effects of the geometry could be important at scales where the matter density was not nec- 
essarily high. In that case, in the semiclassical regime the physical results deviated significantly 
from the predictions made by general relativity [20]. Improved dynamics solves this problem. 
Furthermore, it has been proved that it is the only minisuperspace quantization (among a cer- 
tain family of possibilities) yielding to a physically admissible model [92], independent of the 
fiducial structures, with a well-defined classical limit in agreement with GR, and giving rise to 
a scale of Planck order where quantum effects are important and solve the singularity problem. 



3.3.2 Representation of the Hamiltonian constraint 

When trying to promote the gravitational part of the scalar constraint (3.1) to an operator, we 
find an additional difficulty concerning the inverse of the volume, 

1 _ ^f°q 

V ~ ^/\det{E)\V ' 

The volume operator has a discrete spectrum with the eigenvalue zero included, so its inverse 
(obtained by using the spectral theorem) is not well defined in zero. Nonetheless, following LQG 
[185, 187], from the classical identity 

f fJ=S =E^^7^°^trK(c){[/ l [( c) r 1 ,y}r l ), (3.6) 

we can obtain an operator for the left-hand side of this expression by promoting the functions 
on the right-hand side to the corresponding operators, and by making the replacement { , } — > 
— (i/h) [" , "]■ Note that the parameter I labels a quantization ambiguity. In order not to introduce 
new scales in the theory, we take for I the value jl = ^ A/\p\ [21]. 

Plugging this result into the Hamiltonian constraint (3.1), as well as the curvature given in 
equation (3.4), we obtain that the geometry (or gravitational) contribution to the Hamiltonian 
constraint operator is [21] 

4rav = i^^^V [sml^)sgn^)] 2 (AfpYAf-p, - N-pVMp) , (3.7) 



with 

sin(/2c) 



N 2f i - N- 2fi 
2i 



Let us now deal with the representation of the matter contribution, given in the second term 
of equation (3.2). To represent the inverse of the volume, we follow the same strategy as before, 
now starting with the classical identity 

sgn(p) 1 1 



As before, we take the trace in the fundamental representation and we choose I equal to fl in the 
quantum theory. To fix the ambiguity in the constant s > 0, we choose for simplicity s = 1/2. 
We obtain 
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The action of this operator on the basis states is diagonal and given by 



\v) = b(v)\v), b(v) 



2(27r 7 /2 1 VA) 1 /3 l 



|l/3 I 



1/3 



\v - 1 



1/31 



While, for large values of v, b(v) is well approximated by the classical value l/y/\p\, for small 
values of v they differ considerably. In fact, the above operator is bounded from above and 
annihilates the zero- volume states. 

The matter contribution to the constraint is then given by the operator 





r) 2 






i 


y 




V 




_VW\. 



C mSL t — Sirlpih 



In order for the Hamiltonian constraint operator C = Cg rav + Cmat to be (essentially) self- 
adjoint, we need to symmetrize the gravitational term (3.7). There is an ambiguity in the chosen 
symmetric factor ordering and several possibilities have been studied in the literature [12, 21, 127, 
144, 154, 203] (see [154] for a detailed comparison between them). Due to its suitable properties, 
here we will adopt the prescription called sMMO in [154] 6 , that is a simplified version of the 
prescription of [144]. Its two main features are: 

i) decoupling of the zero- volume state \v = 0); 

ii) decoupling of states with opposite orientation of the densitized triad, namely states \v < 0} 
are decoupled from states \v > 0). 

As we will see, this will give rise to simple superselection sectors with nice properties. Re- 
markably, the behavior of the resulting eigenstates of the gravitational part of the constraint 
already shows the occurrence of a generic quantum bounce dynamically resolving the singularity. 
Therefore, this prescription ensures that the quantum bounce mechanism is an intrinsic feature 
of the theory, independent of the particular physical state considered 7 . 
Then, following [144, 154], we take 



C 



1/2 



6 



h 2 + %-kgPI 





' 1 " 


) 


V 



nV2 



where the operator f2 is defined as 

1 r~r 3 /4 



n 



\P\ 



[mp, - A/Lap) sgu(p) + sgn(p) (J\f 2i j - M- 



2fl) 



\P\ 



r 3/4 



(3.9) 



(3.10) 



The action of sgn(p) on the state \v = 0) can be defined arbitrarily, since the final action of Q 
is independent of that choice, provided that Q\0) = 0. 

Thanks to the splitting of powers of p on the left and on the right, C annihilates the subspace 
of zero- volume states and leaves invariant its orthogonal complement, thus decoupling the zero- 
volume states as desired. We can then remove the state |0) and define the operators acting on 
the geometry sector on the Hilbert space % gra v defined as the Cauchy completion (with respect 
to the discrete measure) of the dense domain 

Cyl s = span{| V ); «eK\{0}}. 



6 The acronym "MMO" refers to the model of [144], by Martm-Benito, Mena Marugan, and Olmedo. 

7 In [21], the quantum bounce was shown just for particular semiclassical states. Then, with the factor ordering 
adopted in [12], it was shown that the quantum bounce is generic, but the result is only obtained for a specific 
superselection sector. The results of [144] are instead completely general. 
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As a consequence, the big bang is resolved already at the kinematical level, in the sense that 
the quantum equivalent of the classical singularity (namely, the eigenstate of vanishing physical 
volume) has been entirely removed from the kinematical Hilbert space (see also [40]). 

In view of the operator (3.9), it is more convenient to work with its densitized version, defi- 
ned as 



C 



-1/2 



" 1 " 




' 1 " 




c 




V 


V 



-1/2 



6 



fl 2 + 8vrGPi 



-fPd 2 , become Dirac observables that commute with the 



since the operators ft 2 and P 2 - 
densitized constraint operator C. Note that, if we had not decoupled the zero- volume states, 

zero would be in the discrete spectrum of [1/V] and the operator [1/V] (obtained via spectral 
theorem) would be ill defined. Nonetheless, in % gra v (with domain Cyl s ) it is well defined. Both 
the densitized and original constraints are equivalent, inasmuch as their solutions are bijectively 
related [144]. 



3.4 Analysis of the Hamiltonian constraint operator 

With the aim of diagonalizing the Hamiltonian constraint operator C, let us characterize the 
spectral properties of the operators entering its definition. As it is well known, the operator 
P? = —h 2 d^ is essentially self-adjoint in its domain <S(M), with double degenerate absolutely con- 
tinuous spectrum, its generalized eigenfunctions of eigenvalue (hu) 2 being the plane waves e^'l"^. 
The gravitational operator Q 2 is more complicated and we analyze it in detail in the following. 

3.4.1 Superselection sectors 

The action of Q 2 on the basis states \v) of the kinematical sector % gra v is 

fl 2 \v) = -f + (v)f + (v + 2)\v + 4) + [fl(v) + f 2 _{v)] \v) - f-(v)f-(v - 2)\v - 4), 

where 

f±(v) = —^V\v ± 2\y/\v\s±(v), s±(v) = sgn(v ±2) + sgn(v), 

so that Q 2 is a difference operator of step four. In addition, note that /_(«)/_(« — 2) = if 
v G (0,4] and f + (v)f+(v + 2) = if v G [—4,0). In consequence, the operator Q 2 only relates 
states \v) with support in a particular semilattice of step four of the form 

Cf = {v = ±(e + 4n), n G N}, eG(0,4]. 

Then, Q 2 is well defined in any of the Hilbert subspaces Hf obtained as the closure of the 
respective domains Cyl^ = lin{|f), v G with respect to the discrete inner product. The 

non-separable kinematical Hilbert space T-L gTa , v can be thus written as a direct sum of separable 

subspaces %g rav = ® £ {Ht © ?C)- 

The action of the Hamiltonian constraint (and that of the physical observables, as we will 
see) preserves the spaces %f (g> L 2 (R,d(f)), which then provide superselection sectors. Therefore, 
we can restrict the analysis to any of them, e.g., to Hf <g> L 2 (M, dcp), for an arbitrary value of 
£6(0,4]. 

The fact that the gravitational part of the Hamiltonian constraint is a difference operator 
is due to the discreteness of the geometry representation, and therefore it is a generic feature 
of the theory. Actually, the different factor orderings analyzed within the improved dynamics 
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prescription (e.g., [12, 21, 144]) display superselection sectors having support in lattices of step 
four. The difference between the superselection sectors considered here [144] and those of [12, 21] 
is that the formers have support contained in a semiaxis of the real line, whereas the support of 
the latters is contained in the whole real line. 

3.4.2 Self-adjointness and spectral properties 

Though the gravitational part of the Hamiltonian constraint operator is not a usual differential 
operator but a difference operator, there exists a rigorous proof showing that it is essentially 
self-adjoint [126]. Here we sketch that proof for the operator that we are considering, £l 2 , but 
indeed the proof can be extended for the different orderings explored in the literature (e.g., 
[12, 21]) 8 . 

In [126] the authors define certain operator H' Aps , 9 which is a difference operator of step 
four, and they show that -f^PS * s unitarily related, through a Fourier transformation, to the 
Hamiltonian of a point particle in a one-dimensional Poschl-Teller potential, which is a well- 
known differential operator. In particular, it is essentially self-adjoint, and then so is H' Aps as 
well. 

In our notation, H' APS ^ s defined on the Hilbert spaces: 

• Hf © U A _ £ , with domain Cyl+ U CylJ_ e , if s ^ 4; 

• H~l © H A © Ho, (Ho being the one-dimensional Hilbert space generated by \v = 0)), with 
domain Cyl| U CyP 4 U lin{|0)}, if e = 4. 

Now, one can show that fl 2 and [4/(3irG)]H' APS (defined on the same Hilbert space) differ in 
a trace class symmetric operator [144, 154]. Then, a theorem by Kato and Rellich [131] ensures 
that Q 2 , defined in the same Hilbert space as H' A p S , is essentially self-adjoint. From this result, 
it is not difficult to prove also that the restriction of VL 2 to Hf (the subspace where we have 
restricted the analysis) is also essentially self-adjoint [144], just by analyzing its deficiency index 
equation [165]. 

On the other hand, it was shown in [126] that the essential and the absolutely continuous 
spectra of the operator H' Aps are both [0, oo). Once again, Kato's perturbation theory [131] 
allows one to extend these results to the operator £l 2 defined in Hf ©%J_ £ . In addition, taking 
into account the symmetry of Q 2 under a flip of sign in v and assuming the independence of the 
spectrum from the label e, we conclude that the essential and absolutely continuous spectra of £l 2 
defined in Hf are [0, oo) as well. Besides, as we will see in next subsection, the (generalized) 
eigenfunctions of ft 2 converge for large v to eigenfunctions of the WDW counterpart of the 
operator. This fact, together with the continuity of the spectrum in geometrodynamics, suffices 
to conclude that the discrete and singular spectra are empty. 

In summary, the operator Q 2 defined on Tif is a positive and essentially self-adjoint operator, 
whose spectrum is absolutely continuous and given by M + . 

3.4.3 Generalized eigenfunctions 

Let us denote by |e^) = ^2 veC + e £ x (v)\v) the generalized eigenstates of Q 2 , corresponding to 
the eigenvalue (in generalized sense) A G [0, oo). The analysis of the eigenvalue equation 
J7 2 |e^) = A|e^) shows that the initial datum e|(e) completely determines the rest of eigen- 
function coefficients e|(e + An), n E N + [144]. Therefore, the spectrum of il 2 , besides being 
positive and absolutely continuous, is also non-degenerate. We choose a basis of states |e^) 

8 Q, 2 is analog to the operator Q of [21]. 

9 The acronym "APS" refers to the model of [21] by Ashtekar, Pawlowski and Singh. 
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normalized to the Dirac delta such that (ef |ef /) = <5(A — A'). This condition fixes the complex 
norm of ef (e). The only remaining freedom in the choice of this initial datum is then its phase, 
that we fix by taking e?(e) positive. The generalized eigenfunctions that form the basis are then 
real, a consequence of the fact that the difference operator has real coefficients. In short, 
the spectral resolution of the identity in the kinematical Hilbert space 7if associated with Q 2 
can be expressed as 

1= f d\\e\)(e\\. 

The behavior of the eigenfunctions e\{e) in the limit v — > oo allows us to understand the 
relation between the quantization of the model within LQC and that of the standard WDW 
theory, where a Schrodinger-like representation is employed in the geometry sector, instead of 
polymeric. Let us study this limit. 

In the WDW theory the analog to the operator J7 2 is simply given by [144] 

2 

2 = ~ [l + 4^ + 4(^) 2 ] , 

^2 

where a = Att^I^. ft is well defined on the Hilbert space L 2 (IR + , dv). Moreover, it is essentially 
self-adjoint, and its spectrum is absolutely continuous with double degeneracy. The generalized 
eigenfunctions corresponding to the eigenvalue A G [0, oo) will be labeled with uj = ±vA G M. 
and are given by 

1 / In If I \ 
= /o i i exp ~ IUJ • ( 3 - U ) 

These eigenfunctions provide an orthogonal basis (in a generalized sense) for L 2 (R + ,dv), with 
normalization (e^e^i) = 5{u — to'). 

Using the results of [128], one can show that the loop basis eigenfunctions e £ x (v) converge for 

^2 

large v to an eigenfunction of the WDW analog Q, . The WDW limit is explicitly given by [144] 

e \( v ) r { ex P [i<l>e{w)] e w (f) + exp [-i<t> e ((Jj\ e^^)}, 
where r is a normalization factor. In turn, the phase 4> £ (uj) behaves as [128, 154] 

Mu) = T(\lu\) + c £ + R £ (\u\), 
where T is a certain function of c £ is a constant, and lim i? e (|w|) = 0. 

3.5 Physical structure 
3.5.1 Physical Hilbert space 

We are now in a position to complete the quantization of the model. In order to do that, we 
can follow two alternative strategies: 

• We can apply the group averaging procedure [18, 137, 138, 139, 140]. The physical states 
are the states invariant under the action of the group generated by the self-adjoint extension 
of the constraint operator, and we can obtain them by averaging over that group. In 
addition, this averaging determines a natural inner product that endows the physical 
states with a Hilbert structure. 
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• We can solve the constraint in the space (Cylg ® <S(M))*, dual to the domain of definition 
of the Hamiltonian constraint operator 10 . Namely, we can look for the elements (ip\ £ 
(Cylg ® 5(R))* that verify (tp\C^ = 0. Then, in order to endow them with a Hilbert space 
structure, we can impose self-adjointness in a complete set of observables. This determines 
the physical inner product [166, 167]. 

Both methods give the same result (up to unitary equivalence): the physical solutions are 
given by 11 

■$(v,(f))= / d\e\{v)[^ + {\)e iv ^ + ^-{\)e- iv ^ <t '], (3.12) 
Jo 

where 

" <a) := v3?- 

In addition, the physical inner product is 

/>oo 

(*i|*2> P hys= / dA[$ + (A)^(A)+#_(A)&_(A)]. 
Jo 

Therefore, the physical Hilbert space, where the spectral profiles ip±(X) live, is 
^phys = L 2 (M + , dX) . 



3.5.2 Evolution picture and physical observables 



In any gravitational system, as the one considered here, the Hamiltonian is a linear combination 
of constraints, and thus it vanishes. In other words, the time coordinate of the metric is not 
a physical time, and provides a notion of "frozen" evolution, unlike what happens in theories, 
such as usual QFT, in which the metric is a static background structure. With the aim of 
interpreting the results in a time evolution picture, we need to define what this concept of 
evolution is. To do that, we choose a suitable variable or a function of the phase space, and 
regard it as internal time [132]. 

In the model that we are describing, it is natural to choose <f> as the physical time. In this 
way, we can regard the Hamiltonian constraint as an evolution equation (p. In turn, v plays 
the role of frequency associated to that time. As we see in equation (3.12), the solutions to the 
constraint can be decomposed in positive and negative frequency components 



dXeKv^Xy^t 



that, moreover, are determined by the initial data *&±(v,(fto) via the unitary evolution 



*±(u, (f>) = U±(4>- <h)*±(v, (p ). 
U±U-<fK 



exp 



±t 4 



47r/f >1 /rf 



(3.13a) 
(3.13b) 



10 We do not expect the solutions of the constraint to live in the kinematical Hilbert space Hg <g) L 2 (I 
which is quite restricted, but rather in the larger space 



(Cyl s ®5(E))* 3 HZ ® L 2 (R,d0) 3 Cyl s ®5(R). 



11 See, e.g., [21] for the application of the group averaging method, or [144] as an example of the second method. 
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This allows us to define Dirac observables "in evolution", namely relational observables [101, 
102, 169], and in turn, to interpret the physical results. Let us note first that, in the classical 
theory, although v is not a constant of motion, v(<fr) turns out to be a single- valued function of (fr 
in each dynamical trajectory [21], and then i>|</, =9 j is a well-defined observable for each fixed 
value (fro. It measures the volume at time (fro. The quantum analogue of that observable is the 
operator 

d\^(v,(fr) = l7 + (0-^o)«* + (v,^o) + E7_(0-0o)i;*_(t;,^o). 
We see that, given a physical solution ty(v,(fr), the action of this operator consists in: 

i) decomposing the solution in its positive and negative frequency components, 

ii) freezing them at the initial time (fr = (fro, 
hi) multiplying its initial datum by v, and 
iv) evolving through equation (3.13). 

The result is again a physical solution, and then the operator constructed in this way is 
indeed a Dirac observable. 

Then, the constant of motion = —ihd^ and the operator form a complete set of Dirac 
(and then physical) observables. Note that both the physical observables and the physical inner 
product preserve not only the superselection sectors, but also the subspaces of positive and 
negative frequency. Therefore, any of these subspaces provide an irreducible representation of 
the observables algebra, and the analysis can be restricted, for instance, to the positive frequency 
sector. 

The operator v\ $ allows to analyze the physical results in evolution. Namely, one can compute 
the expectation value of that observable on physical states at different times. We will carry out 
that analysis in the next section for semiclassical states, and see graphically the occurrence of 
the quantum bounce. 

3.6 Dynamical singularity resolution: quantum bounce 

In the classical theory, when the volume of the universe vanishes, the energy density diverges, 
leading to a big bang singularity. Now, in the quantum theory, the structure of the superselection 
sectors, and more specifically the form of the eigenfunctions e\ (v), ensures that the classical big 
bang singularity is replaced by a quantum bounce. Actually, this result is a consequence of the 
following properties: 

• Exact standing- wave behavior: As we have seen, the eigenfunctions e £ x (v) converge in 
the large v limit to a combination of two eigenfunctions of the WDW theory. These 
eigenfunctions, given in equation (3.11), contract and expand in v, respectively, and can 
be interpreted as incoming and outgoing waves. These components contribute with the 
same amplitude to the limit, and in this sense the limit is an exact standing-wave. 

• No-boundary description: On the other hand, the eigenfunctions e\{v) have support in 
a single semiaxis that, moreover, does not contain the putative singularity v = 0. This 
feature is due just to the functional properties of the gravitational operator f2 2 , and not 
derived from imposing any particular boundary condition. In that sense, the eigenfunctions 
verify a no-boundary description 12 . 



12 In quantum cosmology, the concept of no-boundary has been employed in a different sense of the one discussed 
here [112, 113, 114]. 



Introduction to Loop Quantum Cosmology 



19 



These features imply that the incoming component must evolve into the outgoing one, and vice 
versa, since the flux cannot escape through v = 0. Therefore, in the physical solution (3.12), 
restricted for instance to the positive frequency sector, the expanding and contracting compo- 
nents must lead to two branches of a universe, one in expansion and one in contraction, that 
meet at some positive expectation value of v\tj) Q forming a quantum bounce. This result is then 
independent of the considered physical profile ip + (\). 

For semiclassical states in the region of large v, the expectation value of £>|^ is peaked on 
trajectories that show the replacement of the classical big bang by a big bounce, as depicted in 
Fig. 1. This example corresponds to a physical profile given by a logarithmic normal distribution 
of the type 

^+(A) = L^ e -[m(A/A o) ]V(4 CT 2 ) (3 . 14) 

(as the ones considered in [154]), where the parameters A and a are related to the expectation 
value of P<f) and to its dispersion AP^ by the relations [154] 



Around the bounce point, the expectation values approach the classical value very fast, so that 
the semiclassical limit of the quantum theory agrees with general relativity, as desired. Fur- 
thermore, it has been proven for quite a general class of states that semiclassicality is preserved 
through the bounce [93, 128]. 

Another analytic result holds independently of the choice of state, and can be illustrated in the 
b representation [12]. To this purpose, we choose already classically the densitized Hamiltonian 
constraint, i.e., the total Hamiltonian such that the lapse function is equal to the volume, 
N = v / |det(S)| = V. Thus, one avoids the need to rewrite inverse powers of the volume in 
terms of Poisson brackets: 



C = ~E?EH e l \F* + VC mat = -24(vrG) 2 lim ^v 2 sin 2 2b + VC mat . (3.15) 

'y 3 k on A n ^o An 



At the quantum level, one should choose an operator ordering for the Hamiltonian costraint. 
Different orderings correspond to inequivalent definitions of the theory but they may lead any- 
way to very similar physics. In the absence of a guiding principle selecting one particular ordering 
over the others, one can make a choice convenient for calculational purposes. As an example 
illustrating this point, after regularizing equation (3.15) [An — > jl 2 ), choosing the superselection 
sector v = 4n, and quantizing in the b representation (v — > V = —ihdb), the operator ordering 
can be arranged so that 

[3iTll x {s\n2bd b ) 2 - dl] </>] = 0. 

Because one has a discrete one-dimensional lattice in v space and the Fourier transform in b 
space has support on the interval b € (0, 7r/2) [12], one can define 

z := — , In tan b, 

127T/ 2 ! 

so that we get 

{d 2 -d 2 )*[z,4] = 0. 
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Figure 1. Expectation values and dispersions of v\<j, (in red, vertical bars) in the superselection sec- 
tor with e = 1, corresponding to the physical profile given in equation (3.14) with (P^) = 1000 and 
AP l j > /(P f i ) } = 0.1. The quantum evolution is also compared with the classical trajectories, one in expan- 
sion (blue curve, increasing from left to right) and the other in contraction (green curve, decreasing). 
Graph by courtesy of J. Olmedo. 



This expression is formally identical to the Wheeler-DeWitt equation and the ensuing quantiza- 
tion follows step by step [12]. A key difference, however, is that invariance of the wavefunction 
under parity (frame re-orientation) is not gauge-fixed ab initio and physical states are required 
to satisfy ^+[— z, (j)] = — ^f+[z, <f>[- It follows that the left- and right-moving sectors are not 
superselected and must be considered together. In particular, we can write 

<f>] = *+,l[z+] + *+,r[«-] = £(z+) - e(z-), 

where z± = z ± (f> and £ is some function. 

This fact is crucial for the resolution of the big bang singularity. The volume operator in the 
z variable is 

V = —iv* cosh(Koz)d z , 

where v* is a positive constant and kq = y 127rZpj. At any time <p and on any physical state, 
one can show that the expectation value of the volume is 

(\V\) = \V\ |*+) = V* cosh(^), (3.16) 

where > is the minimal volume at the bounce. Equation (3.16) completes the proof that the 
big bang singularity is avoided in minisuperspace LQC. Further evidence comes from noticing 
that matter energy density has an absolute upper bound (approximately equal to 0.41 times the 
Planck density) on the whole physical Hilbert space [12]. We can reach the same quantitative 
conclusion, albeit not as robustly, when looking at the effective dynamics on semiclassical states 
(Section 8). 

3.7 FRW models with curvature or cosmological constant 

In the previous sections, we ignored the contribution both of the intrinsic curvature = (k/2)5* 
and of a cosmological constant A. Here, we sketch scenarios where the universe is not flat 
(k = ±1) and/or A ^ 0. For more details, consult [24]. 



Introduction to Loop Quantum Cosmology 



21 



3.7.1 Closed universe 

The case of a universe with positive-definite spin connection, K = 1, was studied in [22, 70, 
71, 133, 158, 177, 179, 183]. Due to the extra term in the connection, the form of the classical 
Hamiltonian constraint (3.2) as a function of c (related to metric variables as c = ja + K, 
a dot denotes derivative with respect to synchronous time) is modified by the replacement c 2 — > 
c(c— Vo^) + {l+^ 2 )Vo^ In the classical Friedmann equation, this replacement corresponds to 
H 2 —7- H 2 -\-K/a 2 with K = 1, where H := a/ a is the Hubble parameter. The quantum constraint 
and the resulting difference equation are modified accordingly. There is no arbitrariness in the 
fiducial volume V Q , since it can be identified with the total volume of the universe, which is finite 
and well defined. Then, the choice of elementary holonomy is more natural than in the flat case 
and, locally, one can distinguish between the group structure of SU{2) and SO(3) [183]. As in 
the flat case, the constraint operator is essentially self-adjoint [183] and the singularity at v = 
is removed from the quantum evolution [22, 71, 183]. However, instead of a single-bounce event 
one now has a cyclic model [22]. This can be traced back to the fact that the classical and 
quantum scalar constraint have both contracting and expanding branches coexisting in closed- 
universe solutions, while these branches correspond to distinct solutions in the flat case. 

3.7.2 Open universe 

Loop quantum cosmology of an open universe [179, 181, 191] is slightly more delicate to deal with. 
In contrast with the flat and closed cases K = 0, 1, the spin connection is non-diagonal, so that 
also the connection is non-diagonal and it has two (rather than one) dynamical components c(t) 
and C2{t). The Gauss constraint fixes C2 = 1 and one ends up with the same number of degrees of 
freedom as usual. The volume of the universe is infinite as in the flat case, and a fiducial volume 
must be defined. The classical Hamiltonian constraint is equation (3.2) with c 2 — > c 2 — Vo^j 2 - 
The quantum constraint is constructed after defining a suitable holonomy loop; the bounce still 
takes place and the v = big-bang state factors out of the dynamics. 

3.7.3 A/0 

Another generalization is to add a cosmological constant term, positive [71, 129, 158] or nega- 
tive [32, 70, 126, 181]. At the level of the difference equation, these models have been studied 
in relation to the self-adjoint property. 

For A > 0, below a critical value A* (of order of the Planck energy), the Hamiltonian con- 
straint operator admits many self-adjoint extensions, each with a discrete spectrum. Above A*, 
the operator is essentially self-adjoint but there are no physically interesting states in the Hilbert 
space of the model [129]. 

For A < 0, the scalar constraint is essentially self-adjoint and its spectrum is discrete [126] 
(while, we recall, for A = it is continuous and with support on the positive real line), also 
when K = — 1 [181]. As in the A = 0, K = 1 case, the universe undergoes cycles of bounces [32]. 

4 Bianchi I model 

The next step in extending loop quantum cosmology to more general situations consists in 
the consideration of (still homogeneous but) anisotropic cosmologies. The simplest anisotropic 
spacetime is the Bianchi I model, since it has flat spatial sections. This model has been ex- 
tensively studied, owing to its simplicity and applications in cosmology. In fact, prior to the 
development of loop quantum cosmology, its quantization employing Ashtekar variables was al- 
ready analyzed [23, 135, 136]. The first attempts of constructing a kinematical Hilbert space 
and the Hamiltonian constraint operator within a polymeric formalism were done in [40]. Then, 
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soon after the quantization of the flat FRW model was completed within the improved dynamics 
scheme [21], the same programme was applied to Bianchi I, which we shall review now. 



4.1 Classical formulation in Ashtekar— Barbero variables 

For simplicity, we will consider the model in vacuo. Unlike the FRW universe, which is static 
in vacuo, the vacuum Bianchi I model has non-trivial dynamics. Its solutions are of Kasner 
type [130], with two expanding scale factors and the third in contraction, or vice versa. 

Moreover, for later convenience, we will consider a spatial three-torus topology. Therefore, it 
will not be necessary to introduce any fiducial cell, since the model already provides a natural 
finite cell, that of the three-torus, described with angular coordinates {9, a, 5} running from 
to 2vr. 

Like in the isotropic case, we fix the gauge and choose a diagonal flat co-triad °e l a = 5 l a . The 
presence of three different directions requires three variables to describe the Ashtekar-Barbero 
connection and three more for the densitized triad, that is 13 

4 = E? = ^ 5f v^, i = 9,a, 8. 

The Poisson brackets defining the phase space are then {c l ,pj} = 87rG7<5*-. The spacetime metric 
in these variables reads 

4tH \p z pi p\ 
In turn, the phase space is constrained by the Hamiltonian constraint 



2 c e p e c a p a + c e p e c s ps + C a p a C & p8 

C-BI = 2 t7 = U - t 4 - 1 ) 

7 Z V 



In this expression, V = \/\pePaPs\ is the physical volume of the universe. 



4.2 Quantum representation 

In order to polymerically represent this system, we follow the approach described in Sec- 
tion 3.2 [85]. Holonomies hp (c*) = e^ cVl are defined along straight edges of fiducial length 
27r/Xj 6 R and oriented in the fiducial directions, here labeled by i = 9, a, 5. The fluxes of the 
densitized triad through rectangular surfaces of fiducial area and orthogonal to the i-th direc- 
tion, given by E(A l n , f = 1) = [pi/(47r 2 )]j4g, complete the description of the phase space before 
quantization. The configuration algebra is the tensor product of the algebras of quasi-periodic 
functions of the connection for each fiducial direction: Cyl s = ®iCylg = lin{|//0, /v, us)}, where 
the kets denote the quantum states corresponding to the matrix elements of the holonomies 
Affile 1 ) = e l ^ id} ' I 2 in momentum representation. Hence, the kinematical Hilbert space is the 
tensor product % gm v = ®j%g ra v! where %g rav is the Cauchy completion of Cyl s with respect to 
the discrete inner product = Spiv'-' 

The basic operators are pi and M^.. Their action on the basis states is 

Piltn) =Pi(tH)\m), Pi(vi) = ^jipim, A^Imj) = \m + v'i), 

such that [N~pi,Pj] = ih{N^(c l ),pj}. 



In the following, we will not use the Einstein summation convention, unless specified otherwise. 
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4.3 Improved dynamics 

The most involved aspect that one encounters when trying to adapt the quantization of the 
isotropic case to the anisotropic case lies in the implementation of the improved dynamics, 
explained in Section 3.3. In the presence of anisotropies, we need to introduce three minimum 
fiducial lengths /Zj, when defining the curvature tensor in terms of a loop of holonomies. 
Originally, a naive Ansatz was chosen, given by 

i = (4.2) 




This is the simplest generalization of the ansatz of the isotropic case, equation (3.5). As a con- 
sequence, the operators entering the Hamiltonian constraint have the same form as those of the 
flat FRW model. Furthermore, operators corresponding to different directions commute among 
one another. This allows to complete the quantization obtaining the physical Hilbert space [145], 
in the same way as for the FRW model. However, when the topology is non-compact and then 
a finite fiducial cell is introduced, the physical results depend on this fiducial choice [182]. This 
drawback led to the revision of the definition of p,i, and another Ansatz free of these problems 
was proposed, this time given by 14 

(4.3) 

This choice is geometrically better justified (for a discussion about its derivation see [25]). 
Furthermore, this prescription is the only one verifying a remarkable property: For all the 
fiducial directions, the exponents fiiC 1 of the matrix elements M^c 1 ) have a constant and fixed 
(up to a sign) Poisson bracket with the variable 

/ sV\p~eP^Ps\ (A .s 

v = sgn(p e paPs) - — , 9 t= , (4.4) 
27T7/f, 1 vA 

which is proportional to the volume. Note that it coincides with the parameter v (p) of the 
isotropic case if we identify the three fiducial directions. As a consequence, as we will see, the 
volume will suffer constant shifts in the quantum theory, as in the isotropic case. Thanks to 
this property, the improved dynamics prescription (4.3) nicely implements the interplay between 
the anisotropies and the volume. Instead, within the naive prescription given by equation (4.2), 
there is no interplay between the degrees of freedom associated with different directions, because 
of the commutation between the operators acting on different fiducial directions. Thus, apart 
from giving dependencies on fiducial choices, it also seems less physically motivated. 

Because of these reasons, today it is generally accepted that the more correct improved 
dynamics prescription is equation (4.3), which we shall consider in this paper. 



4.4 Hamiltonian constraint operator 

As in the isotropic case, in order to obtain the Hamiltonian constraint operator we cannot 
represent directly its classical form (4.1), but its expression in terms of the curvature tensor. 
For homogeneous models with vanishing spin connection, as Bianchi I, this expression was given 
in equation (3.1). 

For simplicity, we will densitize the Hamiltonian constraint classically, by simply multiplying 
it by the volume V. In this way we avoid the appearance of inverse powers of the volume 

14 Whenever the three indices i, j, k appear in the same expression, we will consider 7^ 0, so that they are 
different. 
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that make the quantum theory complicated. In any case, as seen in the flat FRW model, the 
densitization could be carried out with no problem in the quantum theory. 

In analogy with the isotropic case, but now taking into consideration that the three fiducial 
directions are different, the curvature operator is the quantum counterpart of the classical 
expression 

Fab = "2X> y wl(i J k r' j 6iS k b , h^ h = hfhpihfrHh?)- 1 . (4.5) 

Taking into account equation (4.5), the expression of the densitized triad and the definition 
of p,i, the densitized Hamiltonian constraint for the Bianchi I model reads 

Cbi = -^y 2 ^e^sgn(p i )sgn(p fe )tr(r i ^. fc ). (4.6) 
^ i,j,k 

In order to represent this constraint as an operator, we first need to define the operators J\n. , 
which represent the matrix elements of the holonomies hf l . To define them we follow a similar 
strategy to that adopted in the isotropic case. We start by reparametrizing Pi(m) with a para- 
meter Xiipi) such that the vectorial field ftid^ produces constant translations in Aj. The solution 
is \i(pi) = sgn(pi)y / |pj|/(47r7/p 1 \/A) 1//3 [25]. The difference with respect to the isotropic case is 
that the translations produced by fiid^, being constant with respect to the dependence on Aj, 
do depend on the parameters Xj and Xk associated with the other two directions. In fact, we 
have fiid^ = (2| Xj A fc |) _1 5 Al . 

As in the isotropic case, we define the operator Afp i such that its action on the basis states |Aj) 
is the same as the transformation generated by jiid^ on the parameter Aj, that is 

■A/±/x fl |A0, X a , Xs) = 

and similarly for A/±^ CT and J\±n s • Moreover, inverting the change of variable, we obtain 
Pi\X e ,X a ,X s ) = (4TT 7 ll 1 VA) 2 /hgn(X i )Xf\X e ,X a ,X 5 ). 

As explained in the flat FRW case, one can always choose a suitable factor ordering for the 
Hamiltonian constraint that allows to remove the kernel of the volume operator, generated by 
the states with X$X,jXs = 0. This is what we will consider. Therefore, the operators J\n. are 
well defined. 

The action of M±p, i can be slightly simplified by introducing the variable v defined in equa- 
tion (4.4), that in terms of A's is given by v = 2XgX a X$. Indeed, making the change from, 
e.g., the states \X$, X a , Xs) to the states \v, X a , Xs), one can check that, under the action of J\p i} 
v suffers a constant shift equal to 1 or —1 depending on the orientation of the densitized triad 
coefficients. On the other hand, the variables Ao- and Xs suffer a dilatation or contraction that 
only depends on their own sign and on v (see [25] for the details). 

The variable v is proportional to the volume 

V = \/\P6PaPs\, V\v, X a ,Xg) = 2-Kjll l V~K\v\\v,X (T ,Xs). 

Therefore, as happened in the isotropic case, in this scheme the volume undergoes constant 
translations. The other two variables measure the degree of anisotropy of the system. 

Once we know how to represent the matrix elements of the holonomies, we can promote the 
Hamiltonian constraint to an operator. When symmetrizing it, we will adopt the prescription 



X ff ± 



1 

2|A CT A 5 | 



Xa, Xs 
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of [104], whose factor ordering is analog to that considered in Section 3 (see equation (3.10)). 
Explicitly, it is given by [104, 147] 



Cbi 



1 



(n 6 n a + n a n e + Metis + fis^e + n a n 5 + n 5 ^a) 



(4.7) 



where 



i 



4iVA 



(A/^ -A/l 2 ^)sgn(pi) + sga(p i )(Sf 2 fi i ~ ti-2ik 



V. 



(4- 



This operator differs from that of [25] in the treatment applied to the signs of pi when sym- 
metrizing. Then, Cbi not only decouples the zero-volume states, but also it does not relate 
states with opposite orientation of any of the triad coefficients, namely, states \v, A CT , X§) with 
opposite sign in any of their quantum numbers. Therefore, Cbi leaves invariant all the octants 
in the tridimensional space defined by v, X a and A,5. Hence, we can restrict the study to any of 
them. We will restrict ourselves to the subspace of positive densitized triad coefficients, given by 

CyL} = lm{\v,X a ,X s ); v,X a ,X s > 0}. 

The action of Cbi on the states of Cylg turns out to be 

«) 2 



Cbi\v, A ct , Xs) 



[x-(v)\v-4, X a ,X s )- -x Q (v)\v,X a ,X s )- 



- x£(v)\v, A CT , Xg)+ + x + (v)\v + 4,X a ,X s ) + ]. 
where we have introduced the coefficients 



X-(v) = 2t/v(v - 2)Vv - 4[1 + sgn(f - 4)], x + (v) = x-{v + 4), 
x (v) = 2{v - 2)v[l + sgn(v - 2)], x+(v) = x^ (v + 2), 



(4.9a) 
(4.9b) 



and the following linear combination of states 
\v ± n, Xfj , Xs)± 



v ± n 

v±n, X a , — —Xs ) + 



+ 



v ± n, 
v ± n, 



v±2 
v±2 



v±n 
v±n, — — A ff , 



v 

v±2 



Act , A,5 ) + 



v 



v ± n 
A CT ,— A,) + 



v±2 

4- X V±2 \ 

v±n,X a , A,5 

v 

v ± n 



v ± n, 



v±2 



A(7) 



v±2 



-X 5 



(4.10) 



Note that, in fact, the operator Cbi is we h defined in Cylg , since X-(v) = if v < 4, and 
Xq(v) = if v < 2. Since there is no v = state, the singularity has no longer analog in the 
kinematical Hilbert space, and then is resolved already kinematically. 

4.5 Superselection sectors 

The analysis of the action of Cbi on a generic basis state \v, A*, A£) shows that: 

i) Concerning the variable v, it suffers a constant shift equal to 4 or —4, the latest only if 
v > 4. Therefore, Cbi preserves the subspace of states whose quantum number v belongs 
to any of the semilattices of step four 



£+ = {e + 4fc,fc = 0,l,2,...}, 



£G (0,4]. 
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ii) Concerning the anisotropy variables, A CT and A,5, the effect upon them does not depend on 
the initial quantum numbers A* and Xg, but only on v = e+4fc. Moreover, this dependence 
occurs via fractions whose denominator is two units bigger or smaller than the numerator. 
As a consequence, the iterative action of the constraint operator on \v, A*, Xg), only relates 
this state with states whose quantum numbers A CT and Xg are of the form A a = u> £ A*, 
with u £ belonging to the set 



m n 



e + 2171 s - k " 



„ e + 2n 

m,neN 



k™ € N, z € Z if e > 2, « = if e < 2 k (4.11) 



The set W £ is infinite and, moreover, one can prove that it is dense in the positive real 
line [104]. Nonetheless, it is countable. Therefore, while the variable v has support in simple 
semilattices of constant step, the variables A a take values belonging to complicated sets, but 
they also provide separable subspaces. As a concrete example, we see that, if both e and A* are 
integers, then A a take values in the positive rational numbers. 

In conclusion, the operator Cbi leaves invariant the Hilbert subspaces Hj x * \*, defined as 
the Cauchy completion of 

c <aj,a| = lin {k A <n *s);ve £+, X a = 0J E Xl, UJ £ E W £ , X* € 

with respect to the discrete inner product (v, X a , Xg\v', X' a , Xg) = (W^AaA^AsA' ■ we wm see ) 
physical observables also preserve these separable subspaces Ht\* \*, and therefore they provide 
sectors of superselection, and we can restrict the study to any of them. 



4.6 Physical Hilbert space 

The Hamiltonian constraint operator Cbi is quite complicated and, unlike the isotropic case 
(and unlike previous quantizations of the model [145]), its spectral properties have not been 
determined. Consequently, it has not been diagonalized either. Therefore, the group averaging 
approach is not useful in this situation, and to analyze the physical solutions one has to impose 
the constraint directly on the dual space (Cyl+ A * A *)*. The elements of that space have the 
formal expansion 

M=E Yl Yl ^(v,w £ \*,Q £ \* s )(v,(J £ \*,u s Xt\. 

From the action of Cbi, one obtain that the constraint (^I^bi^ = leads to the following 
recurrence relation, 

ip + (v + 4,X a ,Xs) = — -r^ [xq (v)i/)-(v, X a , Xg) + x^(v)ip+(v, X a , X s ) 
x + (v) 

- x-(v)ij)-{v - 4, X a , Xg)] . 

In this expression, in order to simplify the notation, we have introduced the projections of 
on the linear combinations of six states defined in equation (4.10), namely, 



ip±(v ± n, Xa-, Xg) = (tp\v ± n, X a , Xg)±. 
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Owing to the property a?_(e) = 0, the above recurrence relation, that is of order 2 in the 
variable v, becomes a first-order equation if v = e: 

' l P+(e + 4, X(j, Xs) = — [xq (e)iP— (e, A CT , Xs) + x^(e)^ + (e, A CT , Xs)] . 

Therefore, if we know all the data in the initial section v = e, we obtain all the combinations of 
six terms given by 

i; + (e + 4, A,,, X s ) = V> f e + 4, A,,, + V> f e + 4, ^rp^ A <x, A 5 

+ tp (s + 4, Act, £ -J~ x ^j + tp (e + 4, ^y^A CT , A 5 

. / , e + 2. £ + 4 \ / £ + 4 x e + 2 x 
+ ij) e + 4, Act, — —A* + -0 £ + 4, — — A,,, X s 



£ £ + 2 / \ £+2 £ 

Prom the combinations tp+{e + 4, A^A^), it is possible to determine any of the individual 
terms ip(e + 4, A CT , A^) that compose them, since it has been shown that the system of equations 
that relate the formers with the latters is formally invertible [147]. This has been proven not 
only for v = £ + 4 but also for all v £ Cj~ . In conclusion, the physical solutions of the Hamiltonian 
constraint are completely determined by the set of initial data 

{^{e,X a ,X s )=ip(e,uj £ X*,u) e X* s ), u E ,O e eW £ , A^,A£eIR + }, 

and we can identify solutions with this set. We can also characterize the physical Hilbert space 
as the Hilbert space of the initial data. In order to endow the set of initial data with a Hilbert 
structure, one can take a complete set of observables forming a closed algebra, and impose that 
the quantum counterpart of their complex conjugation relations become adjointness relations 
between operators. This determines a unique (up to unitary equivalence) inner product. 

Before doing that, it is suitable to change the notation. Following [147], let us introduce the 
variables x a = ln(A a ) = ln(A*)+p e . Note that p £ takes values in a dense set of the real line, given 
by the logarithm of the points in the set W e . We will denote that set by Z £ . For each direction 
a = a or 5, we consider the linear span Cyl A * of the states whose support is just one point x a of 
the superselection sector defined by taking the product of A* with all the points in the set Z £ . 
We call T~L\* the Hilbert completion of this vector space with the discrete inner product. 

Then, a set of observables acting on the initial data ip(x a ,xs) := tp(e, x a , xs) is that formed 
by the operators e lXa and U£ a , with p a £ Z £ and a = a, 6, defined as 

e lx °^{x a ,x 5 ) = e ix ^(x a ,x s ), U^4>(x a , x s ) = 4>(x a + p a , x s ), (4.12) 

and similarly for e lXs and U$ s . These operators provide an overcomplete set of observables and 
are unitary in T~L\* ® T~L\* , according with their reality conditions. Therefore, we conclude that 
this Hilbert space is precisely the physical Hilbert space of the vacuum Bianchi I model. 

Owing to the complicated form of the solutions, together with the fact that a basis of eigen- 
functions of the Hamiltonian constraint operator is not known, the evolution picture of this 
model has not been studied yet. It is expected that, as in the isotropic models, a quantum 
bounce solves dynamically the big bang singularity. Indeed, this has been already checked using 
an effective dynamics [84], in the case of the model coupled to a massless scalar field, as in the 
FRW case. In this respect, it is worth commenting that in analyses where a massless scalar field 
is introduced, the latter serves as internal time to describe the notion of evolution, as seen in 
Section 3.5.2. This field (unlike the geometric degrees of freedom) is quantized adopting a stan- 
dard Schrodinger-like representation, which makes straightforward the construction of a family 
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of unitarily related observables parametrized by the internal time. However, in vacuum cases 
such as the Bianchi I model we have just described, where a suitable matter field is not at 
hand, it would be interesting to describe the evolution regarding as internal time one of the 
geometry degrees of freedom. There is a complication because of the fact that such internal 
time is polymeric, since the geometry degrees of freedom are polymerically quantized. Although 
such description has not been carried out for the Bianchi I model within the current improved 
dynamics, it has been nonetheless constructed for the vacuum Bianchi I model quantized within 
the naive improved dynamics given in equation (4.2). That was done in [146], and an analog 
construction could just as well serve to describe the evolution of the current vacuum Bianchi I 
model (quantized within the scheme given in equation (4.3)), using either the volume variable v 
or its momentum as internal time. 

4.7 Loop quantization of other Bianchi models 

Bianchi models are characterized by possessing three spatial Killing fields. In the case of the 
Bianchi I model, the Killing fields commute and, then, it is the simplest of the Bianchi mod- 
els. This is the reason why it has been extensively studied in the literature, and in particular 
in the framework of LQC. Nevertheless, loop quantization has been also extended to other 
Bianchi models (with non-commuting Killing fields), in particular to Bianchi II and Bianchi IX. 
A preliminary loop quantization of these models was already considered just after the birth of 
LQC [40, 59]. More recently, their quantization has been achieved implementing the improved 
dynamics developed for Bianchi I, in [26] for Bianchi II and in [201] for Bianchi IX. In the 
following we summarize the main characteristics of these works. We refer the reader to those 
references for further details. 

Both the Bianchi II and the Bianchi IX models possess non-commuting Killing fields. As 
a consequence, the fiducial triad and co-triad cannot be chosen to be diagonal, a first feature that 
complicates the analysis in comparison with Bianchi I. Moreover, the spin connection of those 
models is non-trivial. After choosing a suitable gauge and appropriate parameterizations for 
the Ashtekar-Barbero variables, the Hamiltonian constraint of both Bianchi II and Bianchi IX 
consists in that of the (analog) Bianchi I model plus an extra term, of course different for 
each model, but that in both cases involves components of the connection and inverses of the 
densitized triad coefficients. 

The kinematical Hilbert space of the Bianchi II and IX models is identical to that of Bianchi I. 
The main difficulty lies in the representation of the components of the connection that appear in 
the extra term of the Hamiltonian constraint 15 . To tackle this issue, the strategy adopted is to 
define an operator representing the connection [26, 202], choosing a suitable loop of holonomies 
and implementing the ideas employed when constructing the curvature operator for Bianchi I. 
As a result, the components c l of the connection turn out to be represented by the polymeric 
operator 

sin \liC l 

Ci = — , 

Mi 

where jii is the minimum length defined in equation (4.3). 

Once the above operator is defined, the representation of the Hamiltonian constraint as 
a symmetric operator follows straightforwardly, in the same way as for the Bianchi I model. 
Once again, classical singularities are avoided in both models, since the kernel of the volume 
operator can be removed from their quantum theories. 



15 The inverse of the components of the densitized triad can be regularized using commutators with holonomies, 
in an analog way as that employed in the FRW model. 
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Part II 

Midisuperspace models in loop quantum 
gravity 

In the previous part we have discussed the loop quantization of several homogeneous cosmological 
models. In all of them the classical cosmological singularity is avoided. Nonetheless, it is natural 
to ask whether the resolution of the singularity is an intrinsic feature of the loop quantization or, 
on the contrary, if it is a result due to the high symmetry of the homogeneous models. In order 
to answer this question, it seems inevitable to extend the loop quantization to inhomogeneous 
systems. Furthermore, it is essential to make this step in order to develop a realistic theory of 
quantum cosmology. Indeed, as the results obtained in modern cosmology indicate, in the early 
universe the inhomogeneities played a fundamental role in the formation of the cosmic structures 
that we observe today. 

The quantization of inhomogeneous models is technically more complicated than that of 
minisuperspaces, since they possess field degrees of freedom, as the full theory. Hopefully, 
facing the loop quantization of midisuperspaces we will get insights about the open problems 
present in loop quantum gravity. 

In this section we will review the status of loop quantization of midisuperspace models. Clas- 
sically, these models have Killing vectors which reduce the degrees of freedom of the metric but 
the number of Killing vectors are low enough to ensure that the remaining degrees of freedom 
are local. The metrics are therefore parametrized by functions of time and spatial coordinate (s). 
Thus, unlike minisuperspace scenarios, midisuperspace models are field theories (for a compre- 
hensive review, see [31]). Questions beyond the reach of minisuperspace models can be addressed 
in the study of midisuperspaces both classically as well as in the context of a particular quan- 
tization scheme. In the context of LQG, these issues may include the construction of Dirac 
observables or of quasilocal observables, the closure of the quantum constraint algebra and, 
most importantly, whether the singularity resolution mechanism of LQC continues to be valid 
in the field theory context. 

There exist a number of midisuperspace models but the loop quantization procedure has been 
attempted only in a few of them so far. One important class of models is obtained by symmetry 
reduction of GR. 

• Spherical symmetry. In 3 + 1 dimensions, a spacetime is called spherically symmetric if its 
isometry group contains a subgroup isomorphic to 50(3), and the orbits of this subgroup 
are 2-spheres such that the induced metric thereon is Riemannian and proportional to the 
unit round metric on S 2 . These are, in a sense, midway between the minisuperspace models 
and models with an infinite number of physical degrees of freedom. Spherically symmetric 
spacetime metrics depend on the radial coordinate, and therefore these models have to 
be treated as field theories. However, in vacuum, the physical solutions are characterized 
by only a single parameter according to Birkhoff's theorem. In that respect they are 
dynamically trivial, although the gauge-fixing procedure is extremely non-trivial. To make 
them into physical field theories, we need to add matter in the form of dust shells. 

Here we are mainly interested in cosmological models while the spherical symmetric models 
are mostly black hole solutions, in particular the Schwarzschild solution. Hence, we will 
briefly mention the progress made in the context of LQG for completeness. 

In [54] the kinematical framework for studying spherically symmetric models in LQG was 
introduced. The volume operator was constructed in [69] and it was shown that the 
volume eigenstates are not eigenstates of the flux operator. Consequently, the standard 
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prescription of constructing the Hamiltonian operator cannot be used. This problem was 
circumvented in [68] , where the Hamiltonian constraint operator was constructed in terms 
of non-standard variables which mix the connection and the extrinsic curvature. This for- 
mulation was extended to explore the question of singularity resolution of Schwarzschild 
black holes in [51] and the Lemaitre-Toleman-Bondi collapse of a spherical inhomoge- 
neous dust cloud in [60] . The choice of variables made in the above programme is similar 
to the polymer quantization of the Gowdy T 3 model, which will be described in detail 
later. Since the basic quantum variables used in these constructions are different from the 
basic quantum variables of the full theory, another approach in loop quantization of the 
spherically symmetric models was explored in [82]. In this approach, the diffeomorphism 
constraint is fixed leaving the Gauss and the Hamiltonian constraints. The latter is then 
applied to the Schwarzschild solution. The exterior solutions agree with the ones obtained 
in geometrodynamics. The interior solution was studied in [83] where it was shown that, 
after a partial gauge fixing, it can be mapped to the minisuperspace Kantowski-Sachs 
model. After loop quantization, the singularity is replaced by a bounce. In [103] the issue 
of the residual diffeomorphism invariance of loop-quantized spherical symmetric models 
has been investigated. Although the quantization programme is still incomplete, the stud- 
ies done so far indicate that the singularity resolution mechanism described in the previous 
section may be a robust feature of loop quantization. 

Another class of midisuperspace models can be roughly classified on the isometry group of 
the metric. In most cases, this is equivalent to a classification based on the number of Killing 
vectors of the metric. 



• Spacetimes with one Killing vector. These models are obtained by symmetry reduction 
corresponding to one-dimensional spatial isometry groups taken to be M or [7(1) . Some 
of the important features of GR are retained, such as diffeomorphism invariance and the 
field-theoretic non-linear nature of the physical degrees of freedom. It was shown in [105] 
that, locally, it is possible to interpret these 3 + 1 dimensional models as 2 + 1 dimensional 
general relativity coupled to two matter fields, a scalar field and a one-form field corre- 
sponding, respectively, to the norm and the twist of the four-dimensional Killing vector 
field. However, not much work has been done so far in loop quantizing this reduced sys- 
tem apart from the preliminary analysis carried out in [121] in terms of complex Ashtekar 
variables. 

• Spacetimes with two Killing vectors. The next level of simplification is to consider two- 
dimensional spacelike isometry groups. Two types of group action and spatial topologies 
have been studied in LQG: 

— Isometry group of M. x U(l) with the spatial topology being M 3 . These correspond 
to Einstein-Rosen cylindrical waves. An preliminary attempt has been made to con- 
struct the kinematic Hilbert space and define the volume operator (see Section 5.11 
of [44]), although there has not been much progress so far. 

— Isometry group ofC/(l)xL r (l) with the spatial topology being T 3 . These correspond 
to Gowdy models. 

If we impose the additional condition that the Killing vectors are mutually orthogonal, 
we obtain polarized models. These are the simplest midisuperspace models with only one 
field theory degree of freedom. An example is the polarized Gowdy T 3 model. This has 
been studied in some detail in LQG from two different perspectives, both of which will be 
described in detail in the next two sections. 
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5 Hybrid quantization of the polarized Gowdy T 3 model 

Gowdy models are among the best known inhomogeneous cosmologies [107, 108]. They repre- 
sent globally hyperbolic vacuum spacetimes, with compact spatial sections and with two spatial 
Killing fields. The simplest example is the linearly polarized Gowdy T 3 model. Indeed, its 
classical solutions are exactly known [125, 161, 162]. They represent gravitational waves prop- 
agating in a closed expanding universe. Its standard quantization was already considered since 
the 70's [33, 34, 35, 91, 94, 159, 164, 189]. Moreover, a Fock quantization of the model, in which 
the dynamics is implemented unitarily, was achieved [88, 89], and it has been shown that this 
quantization is essentially unique [90, 98] (in a sense that we will explain later). 

In order to apply the quantization employed in LQC to this model, the simplest possibility 
is to carry out a hybrid quantization, that combines the polymeric quantization of the degrees 
of freedom that parametrize the homogeneous solutions, with the Fock quantization for the 
inhomogeneities. This hybrid quantization was developed in [104, 141, 151]. Here we summarize 
its construction and main results. 

5.1 Classical description of the Gowdy T 3 model 

The Gowdy T 3 model represents vacuum solutions to the Einstein equations, with two spatial 
Killing fields that commute, and with spatial hypersurfaces homeomorphic to a three-torus. As 
said before, we will consider the linearly polarized model, that possesses an additional symmetry: 
the Killing vectors are orthogonal to hypersurfaces and, therefore, are mutually orthogonal 
everywhere. Letting d a and dg be the Killing vectors, the model admits global coordinates 
{i, 9, a, 5} adapted to the symmetries, with 9, a, 5 G S . 

After a 3+1 decomposition, we can describe the spacetime metric in terms of the three- 
metric q a b induced in the spatial sections foliating the four-dimensional manifold, the densitized 
lapse N = N/^/q, and the shift vector N a , with a, b G {9, a, 5}. Owing to the isometries d a 
and dg, the Gowdy model verifies qg a = = qg$. This condition fixes the gauge freedom 
associated with the momentum constraint in those directions, and implies N a = = N s [152]. 
As a consequence, the metric components only depend on t and 9 and are periodic in the latter. 
This periodicity allows us to decompose the metric components in Fourier modes 16 . On the 
other hand, the condition of linear polarization imposes q a s = 0. Therefore, the three-metric 
is diagonal and can be described by three fields (r, £, 7), that essentially characterize the area 
of the isometry group orbits, the norm of one of the Killing vectors, and the scale factor of the 
metric induced on the set of group orbits. The phase space is then parametrized by those fields 
and by their momenta (P T , P^, P^), and constrained by the ^-momentum constraint and by the 
Hamiltonian constraint [88]. 

In order to prepare the model for its quantization, the gauge is further reduced. One imposes 
that the generator of the conformal transformations, Pj, and the area of the isometry group 
orbits, t, are homogeneous functions. These conditions fix the gauge freedom associated with 
the non-zero Fourier modes of the ^-momentum constraint and of the densitized Hamiltonian 
constraint, and imply that the functions Ng and N are homogeneous [88, 141, 151]. Then, 
two global constraints remain in the model, the spatial average of the ^-momentum constraint, 
generating rigid rotations in the circle, and the spatial average of the densitized Hamiltonian 
constraint. We will denote them by Cg and Cq, respectively. 

16 We adopt the following convention to define the Fourier modes 4> m of a generic field (f)(6): 




m — 




— im0 



32 



K. Banerjee, G. Calcagni and M. Martm-Benito 



The classically reduced phase space can be split into homogeneous and inhomogeneous sectors. 
The homogeneous sector coincides with the phase space of the Bianchi I spacetime with three- 
torus topology. This sector will be quantized a la LQC, and therefore is parametrized by 
the variables {(c\pi),i = 9, a, 5}, with Poisson bracket {c l ,pj} = 87rG7<5*-, as described in 
Section 4. The inhomogeneous sector is given by the non-zero (inhomogeneous) modes of the 
fields unaffected by the gauge fixing, namely { (£ m , Pg m ) , m E Z \ {0}}. This sector will be 
quantized employing the Fock quantization of [88]. To employ this Fock representation, the above 
inhomogeneous modes are in turn described by annihilation and creation variables (a m ,a^), 
defined as those related to a free massless scalar field. This quantization is preferred as long as 
is the only Fock quantization of the deparametrized system in which the dynamics is unitary 
and with a vacuum invariant under the S 1 translations (the remaining symmetry after gauge 
fixing) [90, 98]. 

The three-metric in terms of the chosen variables reads [104] 
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where the inhomogeneities are encoded in the terms 
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On the other hand, the remaining constraints have the following form: 

oo 

Ce = ^2 m(a* m a m - al m a_ m ) = 0, 



m=l 



Cq = Cbi + = 0, Q = G 
In the above expression, 
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and Cbi is the (densitized) Hamiltonian constraint of the Bianchi I model given in equation (4.6). 
In the Hamiltonian constraint, the inhomogeneities appear in the term Hq, that corresponds 
to the Hamiltonian of a free massless scalar field, and in the term i?-L, that represents an 
interaction term. The inhomogeneities are coupled to the homogeneous sector in a non-trivial 
way, so that the feasibility of the hybrid quantization is not straightforward a priori. 
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5.2 Fock quantization of the inhomogeneous sector 

Once the inhomogeneous sector is described with the appropriate annihilation and creation 
like variables, a m and a* m , it is straightforward to get its Fock quantization. With that aim, we 
promote the variables a m and a* m to annihilation and creation operators, a m and (V m respectively, 
such that [a m , fitj = 5 m m- From the vacuum state |0), characterized by the equations 

«m|0) =0, Vm G Z, 

we construct the one-particle Hilbert space, and the associated symmetric Fock space T [197]. 
The annihilation and creation operators are densely defined in the subspace of T given by finite 
linear combinations of n-particle states 

|n) := | . . . ,n_ 2 ,n_i,ni,n 2 , . ..), 

such that s }2 im n m < oo, being n m G N the occupation number (or number of particles) of the 
rn-th mode. We will denote that space by S. Note that the n-particle states provide a basis for 
the Fock space, orthonormal with respect to the inner product (n'|n) — <^n'tv The action of a m 
and aln on these states is 

a m | . . . , n m , • • • ) — \f n m | • • • , n m 1, . . . ) , 
«Ll • • • > n m> • • • ) = Vn m + 1| • • • ,n m + 1, . . . ). 

5.2.1 Generator of translations in the circle 

The constraint that generates translations in the circle, Cg, does not affect the homogeneous 
sector, and then it is represented on the above Fock space. Taking normal ordering, the corre- 
sponding operator is 

oo 

Ce = m ( a ln^m - fil m d_ m ). 

This operator is self-adjoint in the Fock space T . 

The n-particle states annihilated by Cg are those that satisfy the condition 

oo 

^ mX m = 0, X m = n m - n_ m . 

m>0 

They provide a basis for a proper subspace of the Fock space, that we will denote by Ff. 

5.3 Hamiltonian constraint operator 

Physical states must be annihilated as well by the quantum counterpart of the Hamiltonian 
constraint Cg, given in equation (5.2), which involves both homogeneous and inhomogeneous 
sectors. 

In the previous section, we have already described the representation of the inhomogeneous 
sector, with basic operators a m and aln acting on the Fock space J 7 , which thus constitutes the 
inhomogeneous sector of the kinematical Hilbert space. On the other hand, the homogeneous 
sector is quantized following LQC, namely, it is given by the loop quantization of the Bianchi I 
model. As we discussed in Section 4, in the literature two different implementations of the 
improved dynamics has been applied to the Bianchi I model. Therefore, there exist also two 
different descriptions for the hybrid Gowdy model, one adopting the naive Ansatz (4.2) [104, 141, 
151], and another adopting the improved Ansatz (4.3) [104, 147]. Here we will just explain the 
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second description, which adopts the quantization of the Bianchi I model described in Section 4 
when representing the homogeneous sector. This sector of the kinematical Hilbert space will be 
the kinematical Hilbert space H^x* \* > defined in Section 4.5. 

The first term of the Hamiltonian constraint operator, Cq = Cbi + C§, is thus the Bianchi I 
operator (4.7). We just need to construct the operator that couples homogeneous and 
inhomogeneous sector. 

Let us first focus on the inhomogeneous terms. In order to represent the free Hamiltonian Hq 
and the interaction term H-L, defined in equation (5.3), we choose normal ordering. Then, their 
quantum analogs are given by 

oo 

^ mN m , 

N m + %. 



N„ 



m>0 

oo 



int 



E 

m>0 



^ rr, 



m 



both densely defined in the space S of n-particle states. The operator Hq acts diagonally on the 
n-particle states, and then it is well-defined in the Fock space T . On the contrary, H^ nt does not 
leave invariant the domain S. Indeed, the operator Y m annihilates and creates pairs of particles 
in modes with the same wavenumber \m\, and then L creates an infinite number of particles. 
However, one can prove [104] that the norm of HjIJti) is finite for all n G S, and therefore this 
operator, with domain S, is also well defined in the Fock space J-. 

For the homogeneous terms, we recall that the operator fij, defined in equation (4.8), is the 
loop quantum analogue of the classical term c l pi, and that the inverse powers of \pi\ can be 
regularized taking commutators of pi with holonomies. In view of these prescriptions, C% can 
be represented by the symmetric operator [104, 147] 
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The operator Q, so constructed, leaves the sectors of superselection of the Bianchi I model in- 
variant, and then it is in fact well defined on the separable kinematical Hilbert space 'H't \* x +<giJ-. 



5.4 Physical Hilbert space 

In order to impose the Hamiltonian constraint (^(Cg^ = 0, we expand a general state (ip\ in the 
basis of states \v,u e X*,uj £ Xg) of the homogeneous sector. That is, 

(^1=5^ Y Y (v,u} £ X* a ,0 £ X* s \ ® (ip(v,u E \*,u E \$)\, 
where, let us recall, W e is the set (4.11). In the above expression, 



{tp(v,X a ,X s )\ = (^(v,u £ XI,uj £ Xs)\ 
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is the projection of on the state \v,X a ,Xg) = \v,u} £ X*,Co e X^) of the homogeneous sector and, 
in principle, it must belong to the dual space of some appropriate dense domain of the Fock 
space T . 

If we substitute the above expansion in the constraint, and take into account the action of 
the operators affecting the homogeneous sector, we obtain that the projections (tp(v, X a , Xg)\ 
satisfy difference equations in v that, generically, relate data on the section v + 4 with data on 
the sections v and v — 4, as it happened in the Bianchi I model. Following [104], to simplify the 
notation of the resulting equation, we introduce the projections of (ip\ on the linear combinations 
given in equation (4.10). Namely, we define (i/j±(v ± n, X a , X$)\ = (ip\v ± n, A CT , X$)±- Similarly, 
it is convenient to introduce the combinations of states 



\v ±n, X a , X s )'± 



. v ± n . , 
v±n, X a , X s ) + 



+ 



v ± n, 
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v db n 



win v ± 2 
v±n, — — A CT , X s 



v 



-Xa, Xg) + 



v ± n, 



v±2 
v±2 



and to define the projections of on them: (ip±(v±n, X a , Xg)\ = (tp\v±n, X a , Xs)'±. With this 
notation, the solutions of the Hamiltonian constraint satisfy the explicit relation 

ty+(t> + 4, X a , X s )\ - V [be(v, A (j, X s )b* e (v + 4, X a , X 5 )] 2 ^(^ + (v + 4, X a , Xs)\Hf nt 
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where r] 



i \ 2/3 

4t^7aJ ^ s a dimensionless parameter, X a ,Xs) is the function (5.4), and the 
coefficients x±(v) and x (v) were defined in equation (4.9). 

Similarly to the analysis done in the Bianchi I model, it has been investigated whether 
the solution is totally determined (at least formally) by the data in the initial section v = e. 
The presence of the interaction term in the left-hand side of equation (5.5) complicates a direct 
demonstration of the above statement. However, it is possible to obtain such result in terms of an 
asymptotic analysis of the solutions. Note that the model provides a dimensionless parameter rj 
that can be used to develop an asymptotic procedure, without the need to introduce any external 
parameter by hand. This analysis was carried out in [104], and we refer to it for the details. The 
main result of this analysis is that, in fact, the initial data (tp(e, X a , X$)\ (where A CT and A^ run 
over all possible values in their corresponding superselection sectors) completely determine the 
solution. The solutions turn out to be formal, in the sense that the states (ip(v + 4, Ao-, Xg)\ do 
not belong in general to the dual space of 5, owing to the presence of . in their expression. 

The physical Hilbert space can be characterized, even though the solutions are formal. Indeed, 
once we justify that the set of initial data {(tp(e, w £ A*, Co E X^)\ ; uj £ ,ui £ £ W £ } specifies the 
solution, we can identify solutions with their corresponding initial data, and the physical Hilbert 
space with the Hilbert space of such initial data, exactly as we proceeded with the Bianchi I 
model. 

Once again, the reality conditions over a complete set of observables, acting on the initial data, 
univocally determines the inner product that provides the Hilbert structure. Such observables 
are given, for instance, by the overcomplete set of observables of the Bianchi I model, given 
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in equation (4.12), together with a suitable complete set of observables for the inhomogeneous 
sector, given by [104] 

{(a m + ajj ± (a_ m + a)_ m ), i[(a m - ajj ± (a_ m - al m )]; m G N + }. 

These operators represent the real Fourier coefficients of the non-zero modes of the field 
and of its momentum (9), and in fact they are self-adjoint in the Fock space T . 

Finally, imposing the remaining symmetry of translations on S 1 , the result is that the physical 
Hilbert space of the Gowdy model is [104] 

n phys = H K ,x* ® Tf. 

Namely, it is the tensor product of the physical Hilbert space of the Bianchi I model times 
the physical Fock space for the inhomogeneities (defined in Section 5.2). We note that J-f is 
unitarily equivalent to the physical space of the Fock quantization of the deparametrized sys- 
tem [88, 89]. Therefore, the standard quantum field theory for the inhomogeneities is recovered, 
and they can be seen as propagating over a polymerically quantized Bianchi I background. This 
result supports the validity of the hybrid quantization, since this should lead to the standard 
quantization of the system in the limit in which the effects coming from the discreteness of the 
geometry are negligible. This result is not trivial, since the hybrid quantization is introduced 
in the kinematical setting, and the relation between kinematical and physical structures cannot 
be anticipated before the quantization is completed. 

5.4.1 Singularity resolution 

The classical solutions of the linearly polarized Gowdy T 3 model generically display a cosmolog- 
ical singularity [162]. In the parametrization employed for the hybrid quantization of the model, 
this classical singularity corresponds to vanishing values of the coefficients p^. In the quantum 
theory, the kernel of the operators pi is removed and, as a consequence, there is no analog of the 
classical singularity. This resolution of the singularity is kinematical and, therefore, independent 
of the dynamics. It persists in the Hilbert space of the physical states since they do not have 
projection on the kernel of the operators pi. Moreover, they only have support in a sector with 
positive orientation of the coefficients pi and, then, they do not cross the singularity towards 
other branches of the universe corresponding to different orientations. 

A description of the evolution picture of the model is missing, owing to its high complication. 
It is worthy to note that, at least for the choice of the original naive improved dynamics, the 
effective dynamics of the model has been thoroughly analyzed [73, 74]. In particular, it has been 
studied how the inhomogeneities affect the dynamics of the Bianchi I background. Numerical 
simulations show that the effect of the inhomogeneities does not destroy the bounce. For the 
improved dynamics discussed here, a similar analysis has not been done yet, but we can expect 
similar results, since the bounce mechanism appears for both improved schemes. 

5.5 The Gowdy T 3 model coupled to a massless scalar field 

So far we have discussed the hybrid quantization of the linearly polarized Gowdy model in vacuo. 
This model allows, almost straightforwardly, for the introduction of a minimally coupled free 
massless scalar field with the same symmetries as the metric [143]. Indeed, after a suitable 
rescaling of matter modes, these contribute to the constraints Cg and Cq exactly in the same 
manner as the gravitational field £. Also, the Fock quantization of that system (after a complete 
deparametrization) enjoys the same uniqueness results as that of the model in vacuo and, hence, 
there is a preferred Fock description also for the inhomogeneities of the matter field. Therefore, 
the hybrid approach follows exactly in the same way as for the vacuum case. 
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The interest in considering the model filled with matter lies in the fact that FRW-type so- 
lutions are then allowed. Indeed, as we saw before, the vacuum Gowdy model can be seen as 
a vacuum Bianchi I background filled with inhomogeneities propagating in one direction, and the 
subclass of isotropic solutions of the vacuum Bianchi I model represent trivial Minkowski space- 
times rather than flat FRW universes. Nonetheless, in the presence of matter, the flat-FRW 
becomes the isotropic sector of the Bianchi I model. In that sense, there is a subclass of solutions 
of the Gowdy model coupled to matter that can be regarded as a flat FRW background filled with 
inhomogeneities propagating in one direction. Therefore, the linearly polarized Gowdy T 3 model 
coupled to a massless scalar field provides a simple laboratory where to study, at the quantum 
level (by means of the hybrid quantization), interesting physical phenomena such as the back- 
reaction of the (quantum) inhomogeneities on (polymerically quantized) flat FRW cosmologies 
or, vice versa, the effect of the quantum background geometry on the propagation of the in- 
homogeneities [142]. This analysis is intended to be a first step towards a quantum theory 
of FRW plus inhomogeneities. Its character is quite preliminary, since the inhomogeneities of 
Gowdy are just a subclass of the inhomogeneities that one would introduce in the FRW model 
to account for the inhomogeneities that we observe in our universe. Nonetheless, a complete 
quantization of such a system has not been yet achieved, and the hybrid Gowdy model offers 
a suitable setting to start with. Actually, by employing the same hybrid procedure, the flat 
FRW model plus perturbations is being analyzed [153]. The hybrid quantization also applies in 
this more realistic system as long as a unique Fock quantization for the inhomogeneities is at 
hand as well [95, 96, 97]. 

6 Polymer quantization of the polarized Gowdy T 3 model 

In the previous section, we saw a successful quantization scheme of the linearly polarized Gowdy 
T 3 model where the degrees of freedom were split into homogeneous and inhomogeneous sectors. 
The homogeneous sector was quantized using the LQC techniques, while the inhomogeneous part 
was Fock quantized. One of the significant advantages of the hybrid quantization is that the 
calculations are tractable and the tools developed and studied in LQC can be used to address 
questions even in the midisuperspace context. While it is an extremely useful first step in 
quantization of midisuperspace models, it crucially depends on the fact that the inhomogeneous 
degrees of freedom can be treated perturbatively, and it is assumed that there exits a regime 
in which the most important effects emerging from the discretization of the geometry are those 
that affect the homogeneous subsystem. Ideally, we would like to loop quantize the full polarized 
Gowdy T 3 model without separating the degrees of freedom. In this section we review the work 
that has been done so far in that direction, which has been carried out in [28, 29]. 

6.1 Classical theory 

The variables chosen in this section are significantly different from the ones used in the rest of 
the review so far. We shall therefore indicate the steps followed in obtaining these variables. 

6.1.1 Gowdy T 3 model in Ashtekar variables 

In order to loop quantize, we first need to rewrite the Gowdy T 3 model in terms of real Ashtekar 
variables. Canonical quantization of the unpolarized Gowdy T 3 model in terms of the complex 
Ashtekar variables has been given in [122, 150] which we will briefly sketch below in terms of 
the real Ashtekar variables. 

Recall that, owing to global hyperbolicity, spacetime can be decomposed as M. = Et (g> M, 
where is homeomorphic to a three-tours. As in the previous section, let the angular coordi- 
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nates of E| be (9, a, 5), and the two commuting Killing vectors be £f = d a and £f = d$. These 
isometries imply that the Lie derivatives along these two Killing vectors vanish, i.e. 

£& A a = = C&E?. 

The phase-space variables are therefore only functions of 9. The Gauss and the diffeomorphism 
constraint reduce to 

Gi = deE® + eijA 3 a E%, 

V a = {d a A\)Et - {dgAi)El + e) k A{A\E\. 

The vector constraint given by C a = A l a Gi — V a generates spatial diffeomorphisms. 
We now impose the following gauge-fixing conditions: 

E% = = E%, p = <r,8, 1=1,2. 

The constraints Gj and C p are then solved by Ag = = Ap 1 . 

Thus, only one component of the Gauss constraint (G3) and one of the diffeomorphism 
constraint along the 9 direction (Gq =: C) survive together with the Hamiltonian constraint. 
Since none of the quantities depend on a or 5, we can integrate over the torus T 2 and write the 
symplectic structure as 

a 2 r 

n = ^ J d9(dA 3 e A dE e 3 + dA T p A dE p j). (6.1) 

This is the classical phase phase in terms of real Ashtekar variables. One important observation 
is that this is basically a one-dimensional theory. This is useful because in one dimension, un- 
der orientation-preserving coordinate transformations, a tensor density of contravariant rank p, 
covariant rank q and weight w, can be thought of as a scalar density of weight = w + q — p. 
Hence, under a 9 coordinate transformation, E® transforms as a scalar, E^s transform as scalar 
densities of weight 1, A$ transforms as a scalar density of weight 1, and A^s transform as scalars. 



6.1.2 Choice of new variables 

It turns out that these variables are not suitable for loop quantization and we need to make 
canonical transformations similar to those performed for the spherical symmetric case in [54]. 
Note that, for each p, the A T p and Ej rotate among themselves under the U(l) gauge trans- 
formations generated by the Gauss constraint. These suggest that we can perform canonical 
transformations to define the following variables: 

El = E a cos p, E% = E a sin (3, 

E{ = -E s sin p, E s 2 = E s cos p, 

A\ = A a cos(a + P), A\ = A a sin(a + 0) , 

A] = -A s sin(a + jj) , A 2 S = A s cos(a + P) . 

The angles for the connection components are introduced in a particular fashion for later con- 
venience. 

The radial coordinates, E a , E s , A a , As, are gauge invariant and always strictly positive 
(vanishing radial coordinates correspond to a trivial symmetry orbit which is ignored). 



17 In this section, sum over repeated indices is understood. 
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In terms of these variables, the symplectic structure (6.1) gets expressed as 



n 



4vr^ 

«7 



(16 



dAl A dEl + dX A dE a + dY A dE 5 + d/3 A dP? + d/3 A dP^ 



where 

X := Ao- cos(a), y := yl^ cos(a), 

p/3 := sin(o;), P p := -E S A S sin(a). 

It is convenient to make a further canonical transformation: 



p« 



p/3 



pri 



The constraints are greatly simplified and their detailed expressions can be found in [28]. This 
completes the description of the unpolarized Gowdy T 3 model in the variables we have defined. 
The number of canonical field variables is 10 while there is a threefold infinity of first-class 
constraints. There are therefore 2 field degrees of freedom. We now need to impose two second- 
class constraints such that the number of field degrees of freedom are reduced from two to one 
(as it should be in the polarized case). 



6.1.3 Reduction to polarized model 

In terms of the variables defined above, the spatial three-metric is given by 

ds * = cos ^fL_eL^2 + J^L.*^ + J^fL^ 2 - 2-^- sine dadS. 
El cos^E* cos£E s cos£ s 

For the Killing vectors <%- and ds to be orthogonal to each other, the dad5 term in the metric 
should be zero. This implies that the polarization condition is implemented by restricting to 
the £ = sub- manifold of the phase space of the unpolarized model 18 . In order to get a non- 
degenerate symplectic structure, we need one more condition. We expect the two conditions to 
reduce a field degree of freedom. This turns out to be 

x (9):=2P^ + E e 3 d e \n^ « 0. 

Thus, the reduction to the polarized model is obtained by imposing the two polarization con- 
straints 

C~0 5 x «0, {m,x(0')} = 2^5(6-6'). 

We can solve the polarization constraints strongly and use Dirac brackets. Since the polarization 
constraints weakly commute with all the other constraints, the constraint algebra in terms 
of Dirac brackets is same as that in terms of the Poisson brackets and remains unaffected. 
Furthermore, equations of motions for all the variables other than £, P^ also remain unaffected. 
We can thus set the polarization constraints strongly equal to zero in all the expressions and 
continue to use the original Poisson brackets. It also turns out that the basic variables X, Y, 
unlike in the full theory, are actually the extrinsic curvature components in the a and 5 direction, 
respectively. 

18 Actually there are two possible choices, £ = and £ = 7T. We shall take the constraint to be £ = 0, which 
implies £f > 0. 
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The above construction can be carried out equivalently using SU(2) variables by using in- 
dependent r matrices: 



rs{B) 



cos f](9) T\ +sin?7(^)r2, 
— sin r\(9) t\ + co$T](Q)t2, 

t 3 . (6.2) 



The SU(2) formulation is useful in the quantum theory especially while constructing the Hamil- 
tonian constraint operator. 

Let us review the classical phase space we have constructed. For convenience of notation 
we rename E% := £ and Ag := A. The basic configuration variables are X, Y, A, rj and 
the momentum variables are E a , E s , £, P*>, with Poisson brackets of the form {X,E a } = 
(20/77)^5(6 — 9'). 19 The spatial metric is given by 

ds 2 = fLfLde* + s^da 2 + £%rd5 2 . (6.3) 

The expressions of the constraints are greatly simplified: 

G 3 = ^-[de£ + I^], (6.4) 

C e = \- [E a d e X + E s d e Y - Ad e £ + P^dev] , (6.5) 
«'7 



H 



\ (XE a YE 5 + A£(XE U + YE 5 ) + £d e r,(XE a + YE 8 )) - E a T a E 5 T s 



1 1 

k' y/E L7 

i..fae(*e)l 7-Ce 



2k' 1 v/E J 4 ^/E 2 \v/E 



where E = I^I^S 5 . 

It is obvious from these definitions that X, Y, £, r\ are scalars while E u , E s , A, P 71 are scalar 
densities of weight 1, The Gauss constraint shows that A transforms as a £7(1) connection, 
while 77 is translated by the gauge parameter. All other variables are gauge invariant. 

This completes the process of symmetry reduction from the unpolarized to the polarized case. 
This is a consistent symmetry reduction as can be checked by verifying the constraint algebra. 
It is also possible to show that the solutions of the equations of motion are equivalent to the 
standard Gowdy solutions. In this construction, the goal has been to express the polarized 
Gowdy T 3 in terms of variables which are suitable for loop quantization. In particular, they 
allow a simpler choice of edge and point holonomies, a simpler form for the volume operator, 
and also a more tractable expression for the Hamiltonian constraint. Although it may not be 
possible to make the same choice of variables for other midisuperspace models, similar variables 
have been used in the preliminary steps of loop quantization of another midisuperspace model 
which we describe in brief below. 

Plane gravitational waves. The difference between the polarized Gowdy T 3 model and 
the plane polarized (pp) gravitational waves is in the global topology. While the Gowdy model 
has a compact topology, pp waves have the global topology of Minkowski space. The coordinates 
are no longer angular but due to the homogeneity it is possible to choose an arbitrary finite area 
from the plane wavefronts and consider only finite wave packets. The classical phase space in 
Ashtekar variables is constructed in [117, 118] in a similar way as done for the polarized Gowdy 
T 3 model described above. 



3 Only in this section we will absorb the 4ty 2 and use k := 2G/V. 
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Waves travelling only in one direction are considered to avoid the problem of wave collision. 
Finite pulses of pp waves travelling in the positive or in the negative z direction are characterized 
by a null Killing vector satisfying V^fc;,) = 0. This gives rise to two new constraints: 

U+ = E a K a + E 5 K S - d z £, U- = E a K a - E S K 5 - £ In 

E 

It can be also shown that the constraint U- is identically zero on the constraint surface and 
has weakly vanishing Poisson brackets with all the other constraints, i.e. it is gauge invariant 
and conserved under spatial diffeomorphisms and time evolutions. The constraint [/+ Poisson 
commutes with the Gauss, diffeomorphism and Hamiltonian constraints. This constraint can 
be added as a new first-class constraint and the standard constraint algebra can be enlarged. 
This is an additional ingredient in the analysis of pp waves and the system can now be loop 
quantized. 



6.2 Quantum theory 

In this section, we review the loop quantization of the Gowdy model [29]. The methods and 
steps used here closely follow those used in LQG and are to be viewed as first steps towards con- 
structing a quantum theory of the Gowdy model where all the gravitational degrees of freedom 
are loop quantized. 



6.2.1 Basic states 

Since this is a one-dimensional theory, the graphs are just n arcs with n vertices. The configuration 
variable A is a U(l) connection 1-form, so we integrate it along an edge (an arc along S 1 ) and 
by taking its exponential we define the (edge) holonomy variable valued in U(l): 



hi k \A) := exp (i^ Ja\ , keZ. 



The configuration variables and rj E M/Z are scalars and hence no smearing is needed. 

For these we define the point holonomies (at points v) 



hM(X) := exp [i%X(v)\ , h^(Y) := exp [i^Y(v)\ , hfo) := exp [i\ V (v)] , 

where [i, v E R and A E Z. The X, Y point holonomies are interpreted as unitary representations 
of the compact Abelian group Meohri which is the Bohr compactification of the additive group 
of real numbers R. 

The kinematical Hilbert space is thus a tensor product of the Hilbert spaces constructed 
for A, X, Y, T) variables. For A, the Hilbert space can be constructed using U(l) holonomies 
in a procedure similar to full LQG. For X, Y and rj, we can use the point holonomies as in 
minisuperspace LQC, where the quantum configuration space is taken to be the Bohr compacti- 
fication Meohr- By contrast, rj is an angle variable, so the corresponding point holonomy is 
valued in U(l). 

An orthonormal basis on the tensor-product Hilbert space is provided by the "charge network 
functions". They are labelled by a close, oriented graph T with n edges e and n vertices v, 
a U(l) representation k e for each edge, a U(l) representation X v E Z for each vertex, and MBohr 
representations fj, v , v v for each vertex: 

T vwM^ Y ^ II^ (e) ] II »v[h v (X)]v v [h v (Y)]\ v [h v ( V )] (6.7) 

eer vev{r) 
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= n exp ry / ^) II [ exp (*^f x ) exp exp ( iXvT i) ' 

eer ^ Je J veV {Y) 

where V(r) represents the set of vertices belonging to the graph T. Functions where any of the 
labels are different are orthogonal - in particular, two graphs must coincide for non-zero inner 
product. These basis states provide an orthogonal decomposition for the kinematical Hilbert 
space when all the representation labels are non-zero. 

Note that, unlike in the full theory, in this model we have both point and edge holonomies. 
This construction is also significantly different from the hybrid quantization of the previous 
section. There, the loop Hilbert space has only the homogeneous part represented by point 
holonomies similar to LQC, while the inhomogeneous part is Fock quantized. In that case, the 
full Hilbert space is a tensor product of the two. 



6.2.2 Flux operators 

The conjugate variables are represented as 



E a (8) ~ —ijlp\ 



2 5ho(X) d 



where Zpj := k'H. 



5X(6) dh e {Xy 

The flux variables corresponding to E a , E s , P v are defined by integrating these densities 
on an interval X of the circle, eg T a £ := J x E a ,J : 's ) x '■= fz^- &i being a scalar, is already a 
suitable variable. Their actions on the basis functions (6.7) are 

8{9)T V f 7 . , = A fe « +W + fc e -( fl ) r , 



I, 
I 



fpurp _ 7^P1 rp 

' ueV(r)nz 

F S T — V ii T 

r,k,p,,V,\ 2 v T,k,p,,P,\i 

vev(r)ni 
vev(V)ra 



where I is an interval on S 1 . The symbols e =t (0) either refer to the two oriented edges of the 
graph r, meeting at 9 if there is a vertex at 8, or they denote two parts of the same edge if there 
is no vertex at 9. In such a case, the k labels are the same. In case a vertex is an end-point 
of the interval, there is an additional factor of 1/2 for its contribution to the sum 20 . Note that 
classically the triad components E a and E 5 are positive. Fluxes, however, can take both signs 
since they involve integrals which depend on the orientation. We have thus constructed the 
kinematical Hilbert space together with the representation of the basic variables. Next, we turn 
to the construction of composite operators. 



This follows from 



J a 



1 if Xq g (a, b); 
dx 5(x — x ) = ^ 1 if x a = a or xq = b; 

if xq £ [a, b]. 
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6.2.3 Construction of more general operators 

The diffeomorphism covariance requires that all operators of interest are integrals of expressions 
in terms of the basic operators. They also involve products of elementary operators at the same 
point (same 6) and thus need a "regularization" . As in LQG, the general strategy to define such 
operators is: 

1. replace the integral by a Riemann sum using a "cell-decomposition" (or partition) of S ; 

2. for each cell, define a regulated expression choosing suitable ordering of the basic operators, 
and evaluate the action on basis states; 

3. check "cylindrical consistency" of this action so that the (regulated) operator can be 
densely defined on the kinematical Hilbert space via projective limit; 

4. finally, one would like to remove the regulator. 

Since our model is one-dimensional, both the cell-decomposition and the graphs underlying 
the basis states are characterized by finitely many points and the arcs connecting the consecutive 
points. Adapting the techniques used in LQG [14, 186], the products of elementary variables are 
regulated by using a point splitting and then expressing the fields in terms of the appropriate 
holonomies and fluxes. 

A regulator, for each given graph F, then consists of a family of partitions 11^ such that, for 
each e, each vertex of T is contained in exactly one cell. There is also a choice of representation 
labels ko, /jq, uq, Ao which can be taken to be the same for all e. Since each II r can also 
be thought of as being defined by a set of points such that each vertex is flanked by two 
points, any orientation-preserving diffeomorphism will automatically preserve the order of the 
vertices and cell boundaries. Every sufficiently refined partition then automatically becomes 
a diffeomorphism-covariant regulator. We shall assume that the parameter e denoting a family 
of partitions are sufficiently refined and also plays the role of a diffeo-covariant regulator. The 
regulated expressions depend on e and we recover the classical expressions as e — > 0. 

As in LQG, the issue of cylindrical consistency is automatically sorted out by referring to 
the orthogonal decomposition of "Hkini i.e., by specifying the action of the operators on basis 
states with all representation labels being non-zero. A few comments about the subsequent 
construction: 

1. We have assumed the "length of the intervals" to be same and equal to e. This corresponds 
to a "cubic" partition and is chosen for convenience only. 

2. The charges fj, v , v v can take both signs depending on the orientation of the interval. 
However, the eigenvalues of the volume operator must have explicit absolute values. 

3. Xj denotes the i-th cell of the partition. For a given graph, the partition is so chosen that 
each vertex is included in one and only one interval Zj . The intervals which do not contain 
any vertex of the graph, do not contribute to the summation owing to the property of flux 
operators. Hence, the sum collapses to contributions only from the vertices, independent 
of the partition. The action is manifestly independent of e and even though the number 
of intervals go to infinity as e — > 0, the action remains finite and well defined. 

4. Because of this property of the fluxes, we can choose the Q% point in a cell to coincide with 
a vertex of a graph if Xj contains a vertex, or with an arbitrary point if Xj does not contain 
a vertex. 

5. The measure of the integrals in this section is sometimes suppressed for brevity of notation 
but can be clearly understood from the context. 



44 



K. Banerjee, G. Calcagni and M. Martin-Benito 



Volume operator. In the classical expression for the Hamiltonian constraint, powers of 
E := \£\E a E s occur in the same manner as in the full theory. It is therefore natural to consider 
the expression for the volume of a region IxT 2 and construct the corresponding operator. The 
classical volume operator written in terms of basic variables is 



V 



[1 x T 2 ) = / d 3 x^g = 4vr 2 / d6\l\£\E°E & . 

JlxT 2 JX V 



(6.8) 



To obtain the quantum operator, we first rewrite equation (6.8) as a Riemann sum of volume 
of the cells, which we denote as 



1 = 1 J °i 



This has to be written in terms of the flux variables: 



v e (i) *j2 e v\ s ^ E °&)v E5 & 



i=l 



i=l 



d6E° 



dOE 5 



The right-hand side is now expressed in terms of flux variables. The regulated volume operator 
can be defined as: 



v e (i) -.= j2i\m 



4 = 1 



E° 



E* 



x, 



Clearly, this is diagonal in the basis states and its action on a basis state T r g ^ - ^ yields the 
eigenvalue 



Vj 



1 (lll^ 12 



veinv(r) 



\fi v \ \v v \ \k e +( v } + k e -( v ^\ 



(6.9) 



Thanks to our choice of basic variables, the eigenstates of the flux operators are also volume 
eigenstates. 

Gauss constraint. Consider the Gauss constraint (6.4): 



G 3 = [ d9(d e £ + P r ')^ V[/ P , i + £(8 i + e)-£(6 i 
Js 1 ~[ Ux t 

n 



i=l 



PV + £(9i + e)-£(0i) 



Again, this is easily quantized with its action on a basis state T r g ^, 

k e +(v) - k e -(v) 



giving the eigenvalue 



7& E 

veV(r) 



Xv + 



Notice that in the limit of infinitely fine partitions, for a given graph, if there is a vertex v G 2j, 
then there is no vertex in the adjacent cells. As a result, £(6i + i) gives k e +(v)/2 and —£(0j) gives 
—k e ~(v)/2, since 0j divides the same edge and so does 



Introduction to Loop Quantum Cosmology 



45 



Once again, the eigenvalues are manifestly independent of e and the action is diffeo- invariant. 
Imposition of the Gauss constraint can be done simply by restricting to basis states with labels 
satisfying A^ = — (k e +t v ) ~ k e -( v \)/2, Vw E V(T). Since A^ € Z, the difference in the k labels 
at each vertex must be an even integer. We will assume these restrictions on the representation 
labels and from now on deal with gauge-invariant basis states. Explicitly, 



eer L Je J vev(r) 



2 v; 



exp 



We have also used r](v + (e)) — r)(v (e)) = f d$7], where ir^e) denote the tip and tail of the 
edge e. 



6.2.4 Hamiltonian constraint 

Preliminaries. The Hamiltonian constraint is much more complicated. After some manipula- 
tion, we can write (6.6) as a sum of a kinetic and a potential term, 



H : = 
H K :-. 

Hp := 



7 [H K + H P ], 



1 

K 

J d9N(9)^=[XE a YE s + (A + d e r ] )£(XE (7 + YE 6 )]. 



1 



d9N{9) 



[ dON{9) l -d e 
Js 1 2 



- A {d e S) + — [ — 



2£ (d e £) 



d e E i 
E s 



E 



Here it is more convenient to use SU(2)-valued holonomies using the r/-dependent r matrices 
defined equation (6.2): 



~s / d9'A{9') 



he(T) := exp 
K(0) :=exp[/xoX(0)r CT (0)] = cos 
h s (9) := exp [h,Y(0)t*(0)] 



Mo 



X{9) +2r ff (0)sin 



cos 



"0, 



Y{9) +2T 5 (0)sin 



Mo 
L 2 

2 



X(9) 
Y{9) 



.10) 



Each of the 5'C/(2)-valued holonomies, as well as the sine and cosine, are well defined on the 
kinematical Hilbert space. The interval X will typically be a cell of a partition, (9i,9i + e). 
The parameters no, vq are the chosen and fixed representations of MBohn &o = 1 is the fixed 
representation of the U(l). As before the parameter e which denotes a family of partitions, 
also plays the role of a diffeo-covariant regulator. We briefly describe the steps in obtaining 
a well-defined quantum operator; the details can be found in [29]. 

Consider an expression of the form Ti(hihjh^ 1 hJ 1 h} s: {h' k 1 , y/E}), for distinct i, j, k taking 
values 9, <r, 5. For small values of X, Y, f x A, we can make the following approximations: 



K(9){h a (9)-\V(X)} 
h 5 {9){h 5 {9)-\V(l)} 
h e {h d \V(l)}-- 



^7 £(9)J I E S 
' 2 /i ° T<T V(X) 

k' 7 £(9)J x E' t 
^ U0TS V{X) 



k'j E 5 {9)£{9) 
k'j E a (9)£(9) 

^ U0Ts ^eW' 



±1 k Ea k E& 
2 3 V{X) 



k' 7 E u (9)E 5 (9) 
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AfxeA(0), 



E° 



eE°{9), 



E 5 « eE*(0), 



where V(X) is the volume of the interval I. 

Then, the holonomies can be expanded in a power series. Because of the trace, it is enough to 
expand each holonomy up to first order. The surviving terms are quadratic terms arising from 
products of the linear ones and a linear term coming from h^. If one interchanges the i -<->• j 
holonomies, the linear term retains the sign while the quadratic one changes the sign. Thus, 
taking the difference of the two traces leaves us only with the quadratic terms, which are exactly 
of the form needed in Hk- There are derivatives of r\ which arise from the position dependence 
of the r a , Tg matrices: 



Ta{B + e) 

r S (9 + e) 



■ T a {6) « edoT a 
ts{9) « ed e n - 



ed e r)T S (9), 
-edorrr a (9). 



In the quantization of the Hp, we also need to use the following identities repeatedly (in the 
form LHS/RHS = 1): 



Z{I) := e abc Tv [h a {h-\V(l)}h b {h^,V(l)}h c {h~\V(l)}] 

3 



k'7 



fjLou V(T), 

Z a (l) := e abc Tr [h a {h~\ [V {X)) a }h h {h^ , {V {Z)) a }h c {h- 1 , (V(l)) a }] 



3 / K 7 



(6.11) 



.12) 



These are essentially versions of the identity 1 = (| det(e l a )\ / V E) n [188]. 
It is also convenient to define the following families of operators: 



Ol(T,9) := cos 



sin 



6 5 a (Z,9) := cos 



~^X{9) 



;»oX(9) 



' t VoY(9) 



O e a (l,9) 



sin 



cos 



li Y{9) 



V a (l) sin 
V a (X) cos 
V a (I) sin 
V a (l) cos 



m X{9) 



;VoX(9) 



HoY(0) 



HoY(e) 



A )V a (l)sm 



A 



8111 



A) V a {l) cos 



A 



Above, 9 is a point in the interval X and a > is the power of the volume operator. Again, for 
simplicity of notation we will suppress the 9 labels in the above operators. 

The operators := [cos(- • • )V a sin(- • • ) — sin(- • • )V a cos(- • • )], a = 9, a, 5 appear in all the 
terms and are functions of both holonomies and fluxes. To see that this is actually diagonal 
in the charge network basis, write the cos and sin operators as sums and differences of the 
exponentials (i.e., holonomies). It then follows that 

cosf - • • )V a sin(- • • ) - sin(- • • )V a cos(- • • ) = - \ e - i ^Ve +i ^ > - e +*(~)y e -i(- )1 . ( 6 .13) 

It is now obvious that the operators are diagonal and thus commute with all the flux operators. 
Finally, the operator form of Z a (X) can be obtained as 

Z a {T) := e abc Tv{h a [K\V(l) a ]h b [ h~\ V(xf]h c [ K\V(lf}} 
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i2di{z)d 5 a (z)d e a (i). 



Having noted the ingredients common to the quantization of the different pieces of the Hamil- 
tonian constraint, we turn to each one in some detail. 

Quantization of Hk- Choosing a partition of S 1 with a sufficiently large number of n 
points at Oi, i = 1, ... ,n, 9 n = 2tt, e = 9i + \ — 9i, we write the integral as a sum, 



1 



7 2 ^ WE 

i n 1 r 



r 
+ 



XE a YE 5 + {A + dei]) £{XE C + YE 5 ) (0* 

eAy^i I t: 

J 



Xi 



A + d e r ] )£(9 



X(9i 



E a + Y(9i) / E° 



For small values of the extrinsic curvature components (~ X, Y, classical regime) and sufficiently 
refined partition (e< 1, continuum limit), the i-th term in the sum can be written in terms of 
the traces of the SC/(2)-valued holonomies. The expression in terms of holonomies and fluxes 
goes over to the classical expression in the classical regime. It can be promoted to an operator 
by putting hats on the holonomies and fluxes and replacing Poisson brackets by (ifr) -1 times 
the commutators. Here, the standard choice of putting the holonomies on the left is made. 
Then, we use the expressions for the holonomies in terms of the trigonometric operators given 
in equation (6.10), and evaluating the traces we get the quantum operator as 



A 1 / 

-i-^ y>$) {sin [mjX(9i)] sin W)l } x 



i=l 



2 sin 
2 sin 



v Y(9i + e) 
fi X(9i + e) 



cos 



COS 



u Y(9i 



;MX( 



sin 



sin 



A — A 



i, 



A- A 




where Aj := T}(9i) — rj(9i + e) is outside the integrals. 

Quantization of Hp. All the three terms of Hp are functions only of the momenta, 
but there are a couple of obstacles in a straightforward transcription of Hp. These have to 
be expressed in terms of basic variables, i.e., in terms of fluxes and holonomies. Also, the 
power(s) of momenta in the denominators will make the action on some states singular. The 
first part is easy to take care of thanks to the density weight 1. For the second part, we use the 
identities (6.11) and (6.12) 21 . 

The common strategy followed for these terms is: 

1. introduce a sufficiently large number k > of positive powers of 

1 = \^{K , 1 f m ]- 1 Z{Z)/V{Z), 

and express Z in terms of Z a . This introduces further powers of the volume; 

2. choose a(k) such that explicit multiplicative factors of the volume become 1 and further 
choose k. 



21 In this one-dimensional model, this procedure is equivalent to the point-splitting procedure of [188]. 
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Now the expression can be promoted to an operator. 

First term of Hp. We first rewrite this in terms of the basic variables: 



1 



[ d6N{6) 

7S1 yfE[&) L 4 



=1 



m) 4 e° s Xi E& 



We now follow the strategy mentioned for the terms in the denominator. After some manipula- 
tion, the right-hand side becomes 
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The choice a := 2/3 — l/(3k) removes explicit factors of the volume. The choice of k > is 
limited by a > (being it a power of the volume appearing in Z a ). Some convenient choices 
would be k = 1 (a = 1/3), k = 2 (a = 1/2), and so on. For all of them, the above expression 
can be promoted to a well-defined operator. 

Second term of Hp. To begin with, one observes that E s /E a is a scalar and dg ln(E s / E a ) is 
a scalar density. This term is then manipulated as: 
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We write the right-hand side in terms of flux variables as 
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Now we have the fluxes in the denominator which can be defined exactly as the inverse triad 
operators of LQC. Denoting the fluxes as F c ,x := f x E a ,J : 's i x '■= fx^ 6 , 
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and similarly for T^ x . These can be promoted to well-defined operators. Then, following our 
strategy for the inverse volume factors, we get 

RHS = -\ £ n ^Wt \^k T ^ - tt**) 
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The choice of a would be same as that in the first term. 
Third term of Hp. We can rewrite 
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as 
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where the choice of a is as before. 

We have expressed Hp in terms of the holonomy-flux variables. Quantization can be carried 
out simply via the replacement (Z a ) k — > (—i/h) 3k (Z a ) k . 

We have thus managed to write the Hamiltonian constraint as a well-defined operator on the 
kinematic Hilbert space. Obviously, there are operator ordering ambiguities in the quantization 
of Hk. However, because of equation (6.13), there are no ordering ambiguities in Hp. It is 
also straightforward to verify that the action of the operator is well defined on the states in the 
Hilbert space [29]. 



6.2.5 Ambiguities in the quantization scheme 

In the course of the above construction, there have been obvious issues of quantization ambigui- 
ties. These are in fact ambiguities in the ordering of the operators, in the transcription in terms 
of basic quantum variables as well as in the choice of partitions. 
Let us review them and the choices we made. 

1. In Hk we chose to keep the holonomies to the left. Then, the term containing the volume 
operator acts first on the states and the pieces which do not have any vertices give zero. 

2. In the regularization of H^ we used the inverse volume and plaquette holonomies. We 
could have introduced inverse flux operators and £ operators to replace 1/yE and also 
replaced the X, Y, j x A by wx(^X) / )j,o, and similarly for the others. Such a replacement 
would still give the classical expression back in the limit of small X, Y, e. The quantum 
operator, however, would be different. 

3. The second term in the Hp could be manipulated in terms of inverse powers of y/~E instead 
of introducing inverse flux operators (e.g., by replacing 1/E a = £E S /{y/E) 2 ). This would 
lead to £ 2 CFa^^Fs^+i ~ ^ r 5,x l ^ r o-,x i+ i) and a(k) = 2/3 — 5/ (3k). In the limit of infinite 
refinement, each cell would contain at most one vertex and the cells adjacent to such a cell 
would always be empty. Consequently, the second term of Hp, regulated in the above 
manner, would always give a zero action. 

4. Over and above these different transcriptions, we also have the ambiguities introduced 
by the arbitrary positive power k (and a(k)) and by the arbitrary power I S (0, 1) in 
the definition of inverse flux operators, which is similar to the one in the minisuperspace 
models described before. 

There are also issues related to the choice of partitions, the subsequent e — > limit, and 
the presence/absence of local degrees of freedom. This is most dramatically brought out by the 
second term of Hp. Classically, this is the term which reveals spatial correlations in a solu- 
tion spacetime through dg\n(E s / E a ) [28] and reflects the presence of infinitely many physical 
solutions. In the (vacuum) spherically-symmetric case, such a term is absent and so are local 
physical degrees of freedom. We would like to see if there is a quantization of this term which 
reflects these correlations. This can only be ensured if we chose a partition such that every 
cell has exactly one vertex. Then, the contributions will explicitly depend upon /x, v labels of 
adjacent vertices and, in this sense, spatial correlations will survive in the constraint operator. 
However, the price to pay for that is that we cannot take infinite refinement (e — > 0). 
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An even more restrictive choice would be to pick the partition defined by the graph itself - 
cells defined by the edges and the boundary points of cells as vertices. In this case, the new 
vertices created by Hk would be the already present vertices and the constraint equation would 
lead to a (partial) difference equation among the labels. The e — > limit may then be thought to 
be relevant when states have support on graphs with a very large (but finite) number of vertices; 
heuristically, for semiclassical states. However, more work needs to be done to determine the 
validity of this proposal. 

This completes the kinematic framework of the polymer quantization of this model. Due to 
the complicated nature of the expressions, not much progress has been made beyond that so far. 
A few directions in which there is ongoing further research are: 

• verification of the quantum constraint algebra, 

• obtaining the spectrum and checking the self-adjointness of the Hamiltonian constraint 
operator, 

• construction of Dirac and quasi-local energy observables, at least on the kinematic Hilbert 
space, 

• exploring the possibility of /2-type quantizations in the homogeneous directions, 

• checking whether the Bianchi I cosmological model can be viewed as a sector of this model. 

7 Comparison with the hybrid quantization 

As we have seen in the previous two sections, the polarized Gowdy T 3 model has been quantized 
in two different ways within the framework of LQG. While it would be good to have a procedure 
to compare the physical results of the two frameworks, it is not possible currently. Although 
both quantization schemes start from the classically-reduced phase space of the model, there 
are significant differences in the treatment, both in the classical and the quantum theory. 

The motivation of the hybrid quantization programme is to utilize the tools of LQC on midi- 
superspace models. It tries to determine whether the singularity resolution in LQC is a feature 
of the quantization scheme or an artifact of the high degree of symmetry of minisuperspace. 
Moreover, it provides a suitable arena to analyze the back-reaction between inhomogeneities and 
quantum background geometry. On the other hand, the polymer approach tries to construct 
a loop-quantized theory ab initio, trying to mimic the procedures of the full theory. In addition 
to the fate of the classical singularity in polymer quantization, its aim is to provide a toy model 
where some problems of the full theory, such as the verification of the quantum constraint 
algebra and the construction of observables, can be explored. 

Let us review the progress made in the two quantization schemes so far. 

• The hybrid quantization scheme employs the machinery developed in minisuperspace LQC 
described before, in order to study the midisuperspace model. After a partial gauge fixing, 
it can be easily seen that the polarized Gowdy T 3 model can be thought of as a Bianchi I 
model filled with inhomogeneities propagating in one direction. This fact is exploited to 
break up the degrees of freedom into a homogeneous Bianchi I part and an inhomogeneous 
scalar field. In addition, as a consequence of the partial gauge fixing, only two global 
constraints remain in the model. The polymer quantization programme, on the other 
hand, aims to implement the loop quantization programme by defining suitable Ashtekar 
variables for the entire model. Unlike in the hybrid approach, the system is symmetry 
reduced, but no further gauge is fixed at the classical level, so that the constraints are not 
global but depend on the point. 
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The difference between the two approaches can be seen, on the one hand, in the different 
way the spatial metric is parametrized, by comparing equations (5.1) and (6.3) 22 and, on 
the other hand, in the different constraints surviving in the classically-reduced model. 

• The quantization that is subsequently carried out is also different. In hybrid quantiza- 
tion, the homogeneous Bianchi I is loop quantized while the inhomogeneous scalar field is 
Fock quantized using creation/annihilation operators. The full kinematic Hilbert space is 
a tensor product of the two. There is a non-trivial interaction term in the Hamiltonian 
which couples the homogeneous and inhomogeneous modes. However, all the constraints 
can be expressed as densely defined operators on the tensor product Hilbert space. In the 
polymer quantization, all the degrees of freedom are loop quantized but, unlike the full 
theory, there are both point and edge holonomies. The techniques of full LQG are used 
to define the kinematic polymer Hilbert space. Subsequently, more general operators are 
constructed including the Hamiltonian constraint operator and it can be shown that they 
are well defined on the kinematic Hilbert space and do not depend on the regulator. 

Thus, we have two different quantum theories of the same classical system, both of which 
have a well-defined action of the constraint operators. However, these operators are defined 
on very different Hilbert spaces, and the quantum theories may not be unitarily equivalent. 
While the polymer quantization scheme is closer in spirit to LQG, progress has stalled 
beyond this point because of the extremely complicated nature of the expressions. On the 
other hand, significant progress has been made in the hybrid quantization programme. 

• The construction of the physical Hilbert space has not been carried out so far in the 
polymer scheme. On the other hand, in the hybrid approach the physical Hilbert space 
has been constructed, which turns out to be tensor product of the physical Hilbert space 
of the Bianchi I model and the physical Fock space for the inhomogeneities, which are not 
neglected in this model. Rather, it is possible to view this system as some inhomogeneous 
scalar field on a polymer-quantized Bianchi I space. The classical singularity is absent in 
the physical Hilbert space even in the presence of inhomogeneities. 

Part III 

Effective dynamics 

In the last part, we discuss two aspects of loop quantum cosmology which play an important role 
in the connection between theory and phenomenology. Section 8 presents the effective FRW 
dynamics obtained by evaluating the scalar constraint on semiclassical states. The choice of 
parametrization in minisuperspace is determined by rather robust arguments, but the latters 
undergo several modifications in the context of inhomogeneous models; this is the subject of 
Section 9, where the lattice refinement framework is introduced. 

8 Homogeneous effective dynamics 

The exact and numerical discrete dynamics stemming from the quantum Hamiltonian constraint 
provides important information about the singularity resolution in LQC, but it does not yield 
itself to manipulations suitable for the extraction of inflationary dynamics and observables 
in a semiclassical limit. This can be achieved by evaluating the Hamiltonian constraint on 
semiclassical states, resulting in continuous Friedmann equations corrected by quantum terms. 
At this point, standard analysis techniques developed in classical FRW cosmology can be applied 



Note that there is no relation between the scalar field £ in (5.1) and the angle £ in (6.3). 
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to these equations. The study of linear perturbations (which shall not be reviewed here) requires, 
however, some extra effort. 

Effective cosmological equations of motion are derived from the expression of the Hamiltonian 
constraint on a semiclassical state. The latter is typically decomposed into a gravitational and 
matter sector, |\J/ S c) = YIa b ls rav )A <8> |mat) b- In general, geometrical and matter operators do 
not act separately on physical states because solutions to the Hamiltonian constraint already 
incorporate correlations between the two sectors. So operators on such states are in general 
complicated, entangled observables. However, on a semiclassical state geometrical and matter 
operators commute and they can be treated separately. 

Before discussing how semiclassical states determine an effective dynamics, it is convenient 
to generalize the Hamiltonian constraint and introduce some ambiguity parameters which were 
previously kept fixed. This is done in order to accommodate results which will be later obtained 
in an inhomogeneous setting. In this part we set h = 1. 



8.1 Parametrization of the Hamiltonian constraint 

Let k 2 = 87rG. As before, we define a pair of variables 
fie 6 p 



b :-- 



v :- 
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where p., however, is now an arbitrary dimensionless function of the densitized triad: 
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where n £ R (until now it was fixed to n = 1/2) and p* and a* are, respectively, constants of 
dimension [p*] = —2 and [a*] = 0. Then, 

{b,v} = l. 

In a purely homogeneous model, there is no reason in favour of (and, in fact, there are some 
against) taking n ^ 1/2, but for the time being we do not attempt to justify this generalization. 
Another ambiguity parameter q G R can arise when writing down Thiemann's identity 
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Denote with £q = Aq the area of an elementary plaquette. For a flat homogeneous background, 
the classical scalar constraint becomes [77] 
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The gravitational sector is only a function of b and v. The scalar field part (with potential 
U{(p)) only contains volume factors, 
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As before, the quantum constraint is regularized by assuming that holonomy plaquettes cannot 

1/3 

be shrunk indefinitely, replacing the limit £q — > in equation (8.2) with £$ — > V /2. With this 
substitution, the quantum Hamiltonian operator corresponding to equation (8.2) is well defined: 

C = -4sin26isln26 + C mat , (8.3) 
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We continue to use the notation |v) as the eigenstates of v upon which holonomies act as 
translations, 



v\v) = v\v) 



e iv'b\y\ — \y _)_ w '\ 



These states are also eigenstates of A, 

A\v) = A v \v), 
with eigenvalues 
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8.2 Minisuperspace parametrization 

While in LQG the area spectrum is bounded from below by the minimum area A, due to the 
symmetry reduction the same property is not shared by loop quantum cosmology. Nonetheless, 
we have seen that one may draw inspiration from the full theory and assume that the kinematical 
area of any loop inside the comoving volume V Q is bounded by the area gap for the gauge invariant 
states which are likely to be realized in a homogeneous context. This value is (twice) the LQG 
area gap [25], 

A = 4^/3^7/^, (8.4) 
so that 

(ato) 2 > A. (8.5) 

This step is rather speculative inasmuch as it borrows a result of the background-independent 
framework and forces it into the symmetry-reduced model. It is necessary, however, because 
the quantum scalar constraint in minisuperspace would be singular if one maintained the limit 
£q —> 0. Moreover, the semiclassical limit and the Wheeler-DeWitt equation are reproduced 
correctly. 



54 K. Banerjee, G. Calcagni and M. Martin-Benito 



In a general background, the edges of a spin-network state would intersect a given cell 
only once. By symmetry, the edges of a spin-network state in minisuperspace should traverse 
the fiducial cell, rather than intersecting it from one side; hence the factor of two in equa- 
tion (8.4) [25]. Notice, however, that there is no unique way of fixing the value of A, and 
calculations accounting for different details can produce different numerical prefactors. At any 
rate, these differences are not so large as to give qualitatively inequivalent physical effects. 

If the inequality (8.5) is saturated (smallest possible holonomy path), the comoving cell area 
is also the comoving area gap, that is, the smallest non-vanishing eigenvalue of the area operator 
measuring comoving surfaces. In particular, 

V 2/3 ~ p ~ U M J 
1/3 

One has ]x = £q/V if = A and 

n=\, (8.7) 

a choice corresponding to the improved quantization scheme [12, 21, 22]. The set {\v)} becomes 
the eigenstate basis of the volume operator, v oc p 3 / 2 = V. As the Universe expands, the 
comoving area gap shrinks to zero and the geometry is better and better described by classical 
general relativity, while near the big bang quantum effects become important. 

Originally, the variables p and c were used instead of v and b, corresponding to Jx = 1 (n = 0) . 
In this "old quantization scheme", the states \v) = |/u) coincide with the basis eigenstates of 
the momentum operator p, with eigenvalues v oc p [9, 42]. This case leads to severe restrictions 
of the matter sector if the wavefunctions solving the Hamiltonian constraint are required to be 
normalizable and to reproduce the classical limit at large scales [163]. Also for such reason, 
the improved quantization scheme seems to be the most natural and, as we have already seen, 
the most reasonable in a purely homogeneous context. However, later motivations lead us to 
keep n, p* and the other free parameters of the model as general as possible. In this case, p* is 
some physical squared length determined by the theory which may differ from the mass gap A. 
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8.3 Effective equations of motion 

A semiclassical state l^sc) is peaked around some point (v, b) in the classical phase space. One 
can compute the expectation value of the Hamiltonian constraint operator thereon, using an 
appropriate inner product. Accordingly, for the gravitational part of the Hamiltonian opera- 
tor (8.3) we approximate its expectation value as (^ sc |sin/LicAsin/ic| , I' S c) ~ A v sin 2 fie, and we 
may write (e.g., [44, 100, 184]) 
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where the scalar field energy density is 
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We have two correction functions, a and v. The first is 
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where 



a = 4(1 + n), 



<q<l. 
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The matter correction function is 
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where the ambiguity I lies in the range [67, 188] 
- < I < 1. 
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When a = 1 and the matter sector is a massless free scalar field, equation (8.8) is exact [43]. 
In general, however, the evolution of a finitely-spread semiclassical state will produce quan- 
tum fluctuations leading to additional corrections to equation (8.8) [53, 61]. Assuming that 
the semiclassical wave-packet of the Universe does not spread appreciably, we can stick with 
equation (8.8) also in the presence of a nontrivial scalar potential. Then, the matter energy 
density p is given by equation (8.9). 

The Hamilton equation of motion for the densitized triad gives the Hubble parameter 



H = a 



sin 2/xc 
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In the classical limit, c — > 7a and the right-hand side tends to a/ a for small pc. Combining 
equations (8.8) and (8.14), one gets the Friedmann equation 
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The equation of motion of the scalar field is 
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while the equation for leads to the effective Klein-Gordon equation 
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As v > has a maximum at v = 1 and then decreases down to unity for large v, the friction 
term in equation (8.17) changes sign during the evolution of the universe, the first stage being 
of superacceleration. 

Setting a = 1 = v in the equations of motion (8.15) and (8.17), one ignores inverse- volume 
corrections. On the other hand, in the limit s'm2pc — > 2pc one neglects holonomy corrections 
and the second term in equation (8.15) is dropped. 

The left-hand side of equation (8.15) is positive definite and, if p > (a > if n > —1), the 
energy density is bounded from above: 



p < ap*- 



(8.18) 
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When oc a 2 ( 2n ^ varies with time, there is no constant absolute upper bound. This is avoided 
in the improved quantization (8.7), where the critical density is constant: 

3 

P* ~ ~~ 9 9 • 

Y K P* 

For the particular choice p* = A, the critical density is less than half the Planck density, 
S2ir z 'y 3 

where we used the value 7 ~ 0.238 [148] from the computation of the entropy of non-rotating 
black-hole isolated horizons [7, 8, 148]. The numerical prefactor depends on equations (8.4) 
and (8.5) and it could change in a more complete formulation of the model, but not in a way 
leading to qualitative differences. 

If the ambiguity q is set equal to 1, a = 1 and the lower bound (8.18) is the fixed con- 
stant (8.19) [21, 176, 178]. Thus, the avoidance of the big bang singularity in LQC is confirmed 
at the kinematical level (via the spectrum of the inverse volume operator), by the full quantum 
Hamiltonian dynamics and through the effective dynamical equations. The big bang is replaced 
by a bounce at H = 0, where the energy density is about half the Planck energy. 

These results are encouraging but, of course, insufficient to establish a solid solution of the 
big bang issue. First, the minisuperspace quantization is a toy model. Second, even choosing 
the improved quantization scheme, the critical density ap* may be non-constant if q 7^ 1, in 
which case one might think that the neat bounce interpretation is lost. Nevertheless, we have 
seen how loop quantum gravity could resolve the big bang singularity on one hand, and how far 
we are from a complete understanding of the full theory from the other. 

8.4 Inverse-volume corrections in minisuperspace models 

We discuss now the correction functions a and v from the point of view of their asymptotic 
limits, later stressing an interpretational issue. 

On a semiclassical state the eigenvalues of are approximated by the classical vari- 

able v 1 itself. Consistently, the classical limit corresponds to a large-volume approximation 
where d>1, while in the near-Planck regime ("small volumes"; the reason for quotation marks 
will be soon clear) »<1 Since the momentum operator is v = 2p/p,, the total p-dependence 
of v is effectively 

12^3 p« f]p_ 
a A \p* 

"Near the Planck scale" {v <C 1), the correction functions read 
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From the calculation leading to a and v, one can argue that the "natural" choice of the 
ambiguities I and q can be set at the middle of their range: 

l = \, 1 = \- (8-23) 
In the minisuperspace parametrization, the old quantization scheme corresponds to a = 4 and 

<fa = f, ai = 3T = O(10), <Z„ = ~, v\ = 3^ = O(10 5 ), 



4 



while the improved scheme has 
a = 6 

and 

q a = \, a! = 2 3 / 2 3 3 / 4 » 6, 4, = ^, = 2 5 3 5 / 2 » 500. 

In homogeneous models with n = 0, the duration of this regime depends on the spin represen- 
tation of the holonomies, small j implying a very short super-inflationary period and, actually, 
almost no intermediate stage between the discrete quantum regime and the continuum classical 
limit [67]. Since small- j representations are theoretically favoured, this constitutes a problem. 
It will be relaxed in a different parametrization when inhomogeneities are taken into account. 
In the quasi-classical limit (large volumes), equations (8.10) and (8.12) can be approximated 

as 

a « 1 + a Q 5 Ph (8.24) 
v « 1 + v q 5pi, (8.25) 



where 
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For the natural choice (8.23), the old and improved quantization schemes in minisuperspace 
parametrization correspond, respectively, to 



and 



(j = 4, a = ^ ~ 0.002, vq = $ » 0.02, 



d = 6, a = m ~ °- 01 ' = ifi ~ 0.03. 



Taking g = 1 instead, one gets a negative «o = —1/648 for a = 4 and ao = for (7 = 6. 

Although one can resort to different quantization schemes, equations (8.24), (8.20), (8.25) 
and (8.21) maintain the same structure, where the coefficients a, q a , and q v are robust in the 
choice of the parameters, inasmuch as their order of magnitude does not change appreciably [66]. 
All these parameters can be set to their "natural" values, which are dictated by the form of the 
Hamiltonian or other considerations. 

Now we examine an interpretational issue, already mentioned earlier, related to any paramet- 
rization in pure minisuperspace. On an ideal FRW background, open and flat universes have 
infinite spatial volume and the super-Hamiltonian constraint is formally ill defined because it 
entails a divergent integration of a spatially constant quantity over a comoving spatial slice S, 
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To make the integral finite, it is customary to define the constraint on a freely chosen finite 
region of size V = a 3 V , where V is the corresponding comoving volume: 

/ d 3 x -> / d 3 x = V < +00. 
is Js(v ) 

The volume appears in the correction function (8.22) as 5p\ ~ a~ a ~ V- a / 3 . To make <5 P i 
adimensional, one can use the Planck length lp\ to write 

Physically, the parameter a is related to how the number of plaquettes of an underlying discrete 
state changes with respect to the volume as the universe expands. The latter is a phenomeno- 
logical prescription for the area of holonomy plaquettes, but ideally it should be an input from 
the full theory [45]. For phenomenology at the current level of precision, the most significant 
parameter among {ao,vo,a} is a, which is not as much affected by different choices of the 
minisuperspace scheme. 

Since <5pi is ^-dependent, inverse- volume corrections cannot strictly be made sense of in 
a pure minisuperspace treatment. To cast the problem in other words, the conformal invariance 
of the scale factor a in a non-closed universe make statements such as a* /a <C 1 independent 
of any physical length scale. One could interpret V a as a regulator and send V Q — > 00 at the 
end of calculations, so that in the quasi-classical limit there are no inverse-volume corrections 
at all. However, the full theory does contain these corrections, and one should explain why they 
do not appear in a cosmological setting. At best, this highlights some tension in the theoretical 
construction of the homogeneous LQC effective dynamics. To get a clearer picture, we should 
include inhomogeneities already at the fundamental level. The study of midisuperspace models, 
mentioned in the previous part, is a step in that direction. 

8.5 Models with k ^ and A ^ 

The flat effective dynamical model has been extended to cases with curvature and a cosmological 
constant. 

For a closed universe, K = 1, there is no fiducial volume problem, as mentioned in Section 3.7, 
and inverse- volume corrections are meaningful also in a pure homogeneous and isotropic setting. 
The cyclic bounces appearing in the dynamics of the difference evolution equation [22] exist also 
at the effective level [70, 133, 158]; in particular, the big crunch of classical closed universes can 
be avoided [177]. The bounce persists in an open universe, K = — 1 [191]. In general, all past 
and future strong curvature singularities are resolved in K = ±1 isotropic models; for the closed 
model, weak singularities in the past evolution may also be resolved [179]. 

There is evidence that a cosmological constant, if suitably tuned, does not spoil the singularity 
resolution. When A > and K = 1 [158], the bounce is preserved if the cosmological constant 
is sufficiently small. Above a certain critical value, however, periodic oscillations take place. 
When A < 0, recollapse of the universe is possible, even cyclically [32, 70]. Whatever the sign of 
the cosmological constant, the effective Friedmann equation is equation (8.15), with the critical 
density p* shifted by a constant, A-dependent term. 

9 Inhomogeneous models 

So far we have not given any motivation for taking Ji, oc p~ n . This is the next subject and it 
resides in a framework which does not enjoy the symmetries of a purely FRW background. 
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9.1 Lattice refinement 

In loop quantum gravity, the classical continuum of general relativity is replaced by the ap- 
pearance of discrete spatial structures. It is often expected that the scale of the discreteness is 
determined by the Planck length lp\, but if discreteness is fundamental, its scale must be set 
by the dynamical parameters of some underlying state. Such states are spin networks, graphs 
in an embedding space whose edges e are labeled by spin quantum numbers j e . The quantum 
number determines the area of an elementary plaquette intersecting only one edge e, given by 
A = jlpw/ je(je + !)• The geometrical size of the plaquette changes only when the latter inter- 
sects another edge, thus increasing in quantum jumps. The scale is determined by the Planck 
length for dimensional reasons, but the actual size is given by the spin quantum number. Its 
values in a specific physical situation have to be derived from the LQG dynamical equations, 
a task which remains extremely difficult. However, given the form in which j e appears in the 
dynamical equations, its implications for physics can be understood in certain phenomenological 
situations, such as cosmological scenarios. Then, instead of using the spin labels j e , it is useful 
to refer to an elementary quantum-gravity length scale L, which needs not be exactly the Planck 
length. 

The scale L naturally arises if translation invariance is broken, e.g., by clustering matter or 
inhomogeneous perturbations. The comoving volume V of the system can be discretized as 
a lattice whose J\f cells or patches are nearly isotropic, have characteristic comoving size £q, and 
correspond to the vertices of the spin network associated with V Q . The proper size of a cell is 



To calculate the curvature at the lattice sites within V , we need to specify closed holonomy paths 
around such points. A generic holonomy plaquette is given by the composition of elementary 
holonomies over individual plaquettes. Therefore we set the length of the elementary holonomy 
to be that of the characteristic lattice cell. In other words, the elementary loops of comoving 
size £q we have talked about until now define the cells' walls, while in a pure FRW background 
there is only one cell of volume V Q (the number N is arbitrary). We naturally identify the 
previously ad- hoc function jl(p) as the ratio of the cell-to- lattice size, under the requirement 
that the lattice be refined in time: 



The patch size £q is independent of the size of the fiducial region, since both V D and M scale in 
the same way when the size of the region is changed. Physical predictions should not feature 
the region one chooses unless one is specifically asking region-dependent questions (such as: 
What is the number of vertices in a given volume?). This addresses the issue of conformal 
invariance briefly mentioned above in minisuperspace. In the presence of inhomogeneities there 
is no conformal freedom and, on the other hand, fluxes are determined by the inhomogeneous 
spin-network quantum state of the full theory associated with a given patch [50]. This implies 
that to change the fiducial volume a 3 V would change the number of vertices of the underlying 
physical state. Therefore, there is no scaling ambiguity in the equations of motion [50, 63], 
although the physical observables will depend on the choice of spin-network state. 

The spin-network state described by the lattice can be (and usually is) excited by the action of 
the Hamiltonian operator on the spin vertices, increasing their number and changing their edge 
labels [174, 187]. This process has not yet been established univocally in the full theory, so it is 
convenient to parametrize the number of vertices as in equation (9.1) [45], where the length L(t) 
is state dependent and, by assumption, coordinate independent; its time dependence is inherited 
from the state itself. As the kinematical Hilbert space is usually factorized into gravitational 





/i = AA- 1 / 3 . 
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and matter sectors, the problem here emerges of how to define a natural clock when matter 
does not enter in the definition of a (purely geometrical) spin network. This issue will require 
a much deeper understanding of the theory. So, as unsatisfactory as equation (9.1) may be, we 
take it as a phenomenological ingredient in the present formulation of inhomogeneous LQC. 
The general form (8.1) of p,(p) is obtained if L(t) scales as 



Homogeneous models adopting equation (8.6) feature holonomies which depend on triad vari- 
ables; in other words, curvature components are constrained by the area operator although this 
does not appear in the full constraint. On the other hand, in inhomogeneous models the depen- 
dence of the parameter ft on p is implemented at state (rather than operatorial) level, in closer 
conformity with the full theory [45]. 

As a side remark, the patches of volume L 3 find a most natural classical analogue in inho- 
mogeneous cosmologies, in particular within the separate universe picture [198]. For quantum 
corrections, the regions of size L 3 are provided by an underlying discrete state and thus corre- 
spond to quantum degrees of freedom absent classically. However, the discrete nature of the state 
implies that inhomogeneities are unavoidable and no perfectly homogeneous geometry can exist. 
Given these inhomogeneities and their scale provided by the state, one can reinterpret them in 
a classical context, making use of the separate universe picture. There, the volume V can be 
regarded as a region of the universe where inhomogeneities are non-zero but small. This region 
is coarse grained into smaller regions of volume L 3 , each centered at some point x, wherein the 
universe is FRW and described by a "local" scale factor a(i,x) = a x (i). The difference between 
scale factors separated by the typical perturbation wavelength |x' — x| ~ A <C V 1 / 3 defines a spa- 
tial gradient interpreted as a metric perturbation. In a perfectly homogeneous context, L 3 ~ V 
and there is no sensible notion of cell subdivision of V; this is tantamount to stating that only 
the fiducial volume will enter the quantum corrections and the observables, N = Mq. On the 
other hand, in an inhomogeneous universe the quantity L 3 carries a time dependence which, in 
turn, translates into a momentum dependence. The details of the cell subdivision (number of 
cells per unit volume) are intimately related to the structure of the small perturbations and their 
spectrum. Thus, lattice refinement is better suitable in the cosmological perturbation analysis. 
As long as perturbations are linear and almost scale invariant, the size of the volume within 
which the study is conducted is totally irrelevant. 

9.1.1 Critical density and quantum corrections 

From equations (8.16), (9.2) and (9.1), the critical density is 



In all quantization schemes but the improved one (n = 1/2), is not constant and depends on 
the dynamical patch size L. In any case, the critical density is a number density which depends 
neither on the size of the fiducial volume nor on coordinates, so it is physically well defined even 
outside the improved quantization scheme. 

Similar considerations hold for the quantum correction 8p\. In a purely homogeneous universe, 
the only way to write down equation (8.22) is 5p\ oc (lp\/V 1 ^ 3 ) a , which is volume dependent. 
On the other hand, in the lattice interpretation 



L ~ a' 



3(l-2n) 




(9.3) 



(9.4) 
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and the same quantity is determined by the inhomogeneous state through the patch size L. 
Notice that a > is not the parameter a determined by equation (8.11); n = 1/2 will not 
imply a = 6. The inverse-volume corrections (9.4) do not depend on holonomies due to the 
use of Thiemann's trick (such as equation (3.6)) [57]. Another reason to understand this fact 
is that L 2 is nothing but the expectation value of the flux operator Fg = J s d 2 y Efn a (through 
a surface S with co-normal n a ) on a semiclassical state [57]. In inverse- volume as well as 
holonomy corrections, one refers to elementary building blocks of a discrete state, respectively, 
the plaquette areas and the edge lengths. A pure minisuperspace quantization makes use of 
macroscopic parameters such as the volume of some fiducial region, and fluxes are calculated 

2 /3 

on comoving areas ~ V ■ On the other hand, in the lattice-refinement formulation of loop 
quantum cosmology one uses the microscopic volume of a cell, and fluxes are defined on comoving 
areas ~ £q. This leads to equation (9.4), with some phenomenological parameter a. 

Intuitively, holonomy corrections become large when the Hubble scale H~ 1 = a/ a ~ 7L is 
of the size of the discreteness scale, an extreme regime in cosmology. In terms of the classical 
energy density p = 3H 2 /k 2 , holonomy corrections can be quantified by the parameter 

<W := £ - = (jHL) 2 . 
P* 

These are small when 5hoi "C 1. In order to compare inverse- volume with holonomy corrections, 
we notice that 

5 Pl =( 7 l Pl H5^/ 2 f. (9.5) 

For a universe of causal size H~ 1 ~ lp\, inverse- volume corrections are considerable and be- 
have very differently from what is normally expected for quantum gravity. For small densities, 
holonomy corrections are small, but inverse-volume corrections may still be large because they 
are magnified by an inverse power of <5hoi- As the energy density decreases in an expanding 
universe, holonomy corrections fall to small values, while inverse-volume corrections increase. 
For instance, in an inflationary regime with a typical energy scale of p ~ W~ 10 pp\, we can 
use equation (9.5) with a = 4 to write <5hoi ~ 1CT 9 / 'y/5p\. Small holonomy corrections of size 
<5hoi < 10~ 6 then require inverse-volume correction larger than <5pi > 10~ 6 . This interplay of 
holonomy and inverse- volume corrections can make loop quantum cosmology testable, because it 
leaves only a finite window for consistent parameter values, rather than just providing Planckian 
upper bounds. It also shows that inverse-volume corrections become dominant for sufficiently 
small densities (eventually, of course, they are suppressed as the densities further decrease). 

9.1.2 Lattice parametrization 

The lattice refinement picture allows us to reinterpret minisuperspace quantization schemes 
in a different language. Equation (9.4) replaces the total lattice fiducial volume V as the 
"patch" (i.e., cell) volume L 3 [39]. This means that one makes the formal replacement V — > 
V/N everywhere in minisuperspace expressions, which can be also justified as follows. At the 
kinematical level, internal time is taken at a fixed value but the geometry still varies on the 
whole phase space. In this setting, we must keep N fixed to some constant No while formulating 
the constraint as a composite operator. Since the vertex density does not depend on the choice 
of fiducial volume, it is physically reasonable to expect the Mq factor to be hidden in the 
kinematical quantity a* (or p*). The net result is the Hamiltonian constraint operator of the 
previous sections. 

However, when one solves the constraint or uses it for effective equations, one has to bring 
in the dynamical nature of M from an underlying full state. This is the motivation for promo- 
ting AT to a time-dependent quantity. For some stretches of time, one can choose to use the 



62 



K. Banerjee, G. Calcagni and M. Martin-Benito 



scale factor a as the time variable and represent J\f(a) as a power law (equation (9.2)), 

M = AT a 6n . (9.6) 
Overall, quantum corrections are of the form (9.4), 

S Pi =(ik^Y = (4^)"a^, (9.7) 



where a > 0. This equation cannot be obtained in a pure minisuperspace setting. 

The parameter a plays two roles, one as a dynamical geometric quantity and the other as 
internal time. While writing down the semiclassical Hamiltonian with inverse-volume (and 
holonomy) corrections, one is at a non-dynamical quantum-geometric level. Then, internal 
time is taken at a fixed value but the geometry still varies on the whole phase space. In this 
setting, we must keep N fixed while formulating the constraint as a composite operator. The 
net result is the Hamiltonian constraint operator of the basic formulation of loop quantum 
cosmology [38, 42] not taking into account any refinement, corresponding to n = and a = a. 
On the other hand, equation (9.6) captures operator as well as state properties of the effective 
dynamics. The parametrization of N as a power law of the scale factor is simply a way to encode 
the qualitative (yet robust) phenomenology of the theory. The general viewpoint is similar to 
mean-field approximations which model effects of underlying degrees of freedom by a single, 
physically motivated function. 

Comparing with the earlier minisuperspace parameterization, equation (9.7) gives a = (1 — 
2n)<7. Since dAf/dV > 0, one has n > 0: the number of vertices Af must not decrease with the 
volume, and it is constant for n = 0. Also, £o ~ a 1-2 ™ is the geometry as determined by the 
state; in a discrete geometrical setting, this has a lower non-zero bound which requires n < 1/2. 
In particular, for n = 1/2 we have a constant patch volume as in the improved minisuperspace 
quantization scheme [21]. In contrast with the minisuperspace parametrization (8.11), in the 
effective parametrization of equation (9.7) we have a = for the improved quantization scheme 
n= 1/2. The range of n is then 

< n < -. 

The critical density /?* oc a &<y2n ~ l \ equation (9.3), is still constant for n = 1/2. 

The exponent a in equation (9.6) can be taken as a small positive integer. In fact, the 
correction function <5pi depends on flux values, corresponding to p for the isotropic background. 
Since p changes sign under orientation reversal but the operators are parity invariant, only even 
powers of p can appear, giving a = 4 as the smallest value. Therefore we set a > 4. 

To summarize, a is a time-independent parameter given by the quasi-classical theory and 
with range 



a > 0. 



a may be different in a and v for an inhomogeneous model, but we assume that the background 
equations (8.24) and (8.25) are valid also in the perturbed case. The coefficients «o and vq 
become arbitrary but positive parameters. In fact, from the explicit calculations of inverse- 
volume operators and their spectra in exactly isotropic models and for regular lattice states in 
the presence of inhomogeneities [41, 52, 66], correction functions implementing inverse- volume 
corrections approach the classical value always from above. This implies that 



«o > 0, 



>0. 
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The lattice parametrization replaces the one for homogeneous LQC. In fact, strictly speaking, 
the use of one parametrization instead of the other is not a matter of choice. A perfectly ho- 
mogeneous FRW background is an idealization of reality which, in most applications, turns 
out to be untenable. The study of cosmological perturbations with inverse volume correc- 
tions [55, 56, 57, 62, 63, 64, 65, 81] is an example in this respect. In that case, therefore, the 
lattice refinement parametrization is not only useful, but also required for consistency. Effective 
linearized equations in the presence of holonomy corrections are under development, but we do 
not have complete control over them yet. For vector and tensor modes a class of consistent con- 
straints with a closed algebra is known [62, 63, 156], and therefore inspections of cosmological 
holonomy effects have been analyzed for gravitational waves [86, 109, 110, 155, 157]. On the other 
hand, anomaly cancellation in the scalar sector has been worked out only recently [75, 76, 200]. 

We conclude with some comments on the superinflationary phase of loop quantum cosmo- 
logy. In the near-Planckian regime, the small- j problem in the homogeneous parametrization 
is reinterpreted and relaxed in terms of the lattice embedding. The volume spectrum depends 
on the quadratic Casimir in j representation: p,~ n ~ y 2 / 3 ~ yj Ci(j) ~ j. A higher- j effect 
can be obtained as a refinement of the lattice (smaller fi) [58], thus allowing for long enough 
super acceleration. A change in £i(p) can be achieved by varying the comoving volume V Q . 
This is an arbitrary operation in pure FRW, while in inhomogeneous models fl is a physical 
quantity related to the number of vertices of the underlying reduced spin- network state. As long 
as a calculation of this effect from the full theory is lacking, we will not be able to predict the 
duration of the small- volume regime. More importantly, the closure of the constrained algebra is 
a wide open issue outside the approximation (8.24), (8.25), and early works on superinflationary 
cosmological perturbations [67, 78, 87, 119, 175, 190] have not received a rigorous confirmation. 

10 Conclusions 

The kinematical structure of LQG is well defined although there are some technical difficulties 
in the construction of the physical theory. One way of increasing our understanding of loop 
quantization is to apply it to simple systems. Here we look at applications to cosmological 
systems, also with the hope of making progress in developing a realistic theory of quantum 
cosmology. 

Owing to the presence of an underlying full theory, this procedure has been much more 
successful than the earlier WDW quantization. Trying to follow the steps of LQG, using the 
fact that the underlying geometry is discrete and making physically well-motivated Ansatze 
about the size of the fiducial cells, we have obtained a quantization scheme which ensures the 
resolution of the classical singularity as well as the correct semiclassical limit. This is a generic 
feature of LQC of all minisuperspace models. We have described two examples of minisuperspace 
LQC in detail: 

• The flat FRW model is the simplest and the most rigorously studied model in LQC. We 
have used it to explain the kinematic structure of LQC, and to describe the construction 
of the physical Hilbert space and of observables. These provide an evolution picture of the 
universe with respect to a massless scalar field playing the role of a clock variable, which 
in turn serves to illustrate the mechanism of singularity resolution. 

• Bianchi I model serves as an example for the complications which arise in LQC, from 
the correct choice of the quantization scheme to the construction of the physical Hilbert 
space. In this model the evolution picture is not complete yet as long as a basis of states 
diagonalizing the Hamiltonian constraint remains unknown. The Bianchi I model also 
plays a crucial role in the hybrid quantization of Gowdy T 3 model. 
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Once we have got some handle on the construction of quantum theories of minisuperspace 
models, we need to consider models with field theory degrees of freedom, i.e., midisuperspace 
models. They can serve as good toy models for testing the field-theoretical features of the full 
theory. Also a study of the fate of the classical singularity in these models is needed to ensure 
that the singularity resolution mechanism in LQC is generic and not an artifact of the symmetry 
reduction. However, only one model, the polarized Gowdy T 3 model, has been studied so far 
in LQC. We have described and compared two approaches, the hybrid quantization and the 
polymer quantization procedures. The hybrid quantization scheme, although very successful, 
quantizes a part of the geometry in the Fock way, while the polymer quantization scheme is 
still incomplete. We feel the need to study a complete loop quantization of this and other 
midisuperspace models for a better understanding of LQG/LQC. 

The path we have followed in the first two parts is to bring more and more complicated 
systems under the ambit of loop quantization. We have started with flat FRW, which is ho- 
mogeneous and isotropic. Then we have removed the isotropy condition and studied Bianchi I. 
Finally, we have tried to lift even the homogeneity condition (although retaining some other 
symmetries) and studied the Gowdy T 3 model. 

Instead of looking at various models, one can take another path and look at the generic 
features of the theory and try to incorporate these as corrections to the classical equations. 
This effective equation approach, presented in Part III, is a setting where the effects of LQC 
can be incorporated in a way suitable for phenomenology and observations. We have described 
the progress made in both homogeneous and inhomogeneous contexts. In the latter, the main 
complication is the closure of the effective constraint algebra, whose study has been completed 
for inverse-volume corrections but is still ongoing when holonomy corrections are switched on. 
While verifying this property, we have to ensure that we retain the features which resolve the 
singularity while giving the correct large-scale behaviour. Once this is satisfied, we can use 
this framework to make predictions and constrain the parameters of the model, at least in the 
semiclassical regime. 

A number of issues are under active investigation in Loop quantization of cosmological models. 
Here we mention some interesting problems in the areas of LQC covered in this review. At the 
homogeneous level, in solvable LQC we are still trying to understand the nature of the bounce 
replacing the classical singularity and what are the characteristics of the universe at late times 
after the bounce. This has been extensively explored in the simple models where an evolution 
picture has been developed, but it remains a challenge to control the complete evolution even in 
FRW models with massive scalar fields and in the Bianchi I models. Other Bianchi cosmologies 
need to be analyzed in greater detail because the non-trivial topologies and properties of the 
connections introduce additional subtleties in the quantization. Finding methods for solving 
these issues may indicate ways to address their counterparts in the full theory. Also, it would 
be important to understand whether some of the features of classical evolution (for example, 
the chaotic approach to singularity in Bianchi IX models) are retained after Loop quantization. 
Although we now have a good understanding of the big-bang type singularity from the LQC 
perspective, the nature and the possible resolutions of big-rip type singularities require further 
investigation. 

Here we concentrated on the Hamiltonian formalism and on cosmological models constructed 
via the imposition of quantum constraints. Quantum cosmology can be also defined via different 
sum-over-histories approaches. One is to recast the dynamics of LQC in a form resembling the 
spin-foam formulation of the full theory, and express it in terms of a "vertex expansion" [10, 11, 
80, 116, 120]. This corresponds to a minisuperspace symmetry reduction of the path integral of 
the full theory. Another possibility, also dubbed "spin-foam cosmology" and attempting to get in 
closer contact with the full theory, is to evaluate spin-foam amplitudes on graphs interpreted as 
homogeneous and isotropic geometries [36, 37, 115, 134, 168, 193, 194]. A third line of research, 
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inspired by group field theory, consists in writing down a minisuperspace version of a scalar field 
theory action where wavefunctions are promoted to fields and the Hamiltonian plays the role 
of kinetic operator [79, 106]. Interaction terms, possibly nonlocal, are responsible for topology 
changes and can give effective contributions to the Hamiltonian constraint in a linear mean-field 
approximation. All these approaches are at their first stages of development. 

Much work needs to be done in midisuperspace models where only one cosmological example 
has been studied so far. Although the hybrid quantization scheme has made significant progress, 
one should still fully justify the approach and check whether similar methods are viable for the 
other Gowdy models. Again, that may depend on the development of the other Bianchi setups. 
The polymer quantization scheme provides a genuine toy model suitable for trying to solve the 
problems facing the full theory, and for this reason it should be further refined. We also need 
to look at other midisuperspace scenarios, particularly 1-Killing vector models which can be 
mapped to 2 + 1 gravity with scalar fields. Finally, it is not clear how to obtain minisuperspace 
LQC as a sector of LQG and at what level we expect such an embedding, and whether it is 
possible at all. However, a similar exercise can be attempted in LQC by checking whether one 
can embed minisuperspace models into midisuperspace models. 

The power and importance of the effective-dynamics technique have been slowly gaining mo- 
mentum. Effective dynamics and effective quantum constraints provide the means to ask ques- 
tions such as whether the LQC corrections are consistent with the observed universe and whether 
the underlying discreteness coming from LQC has any appreciable effect on the inflationary 
imprint. All these challenges are currently being undertaken. Confrontation of LQC with ex- 
periments is one of the major open problems to be addressed. Detailed investigations on linear 
cosmological perturbations showed that the size of quantum corrections could be considerably 
larger than what previously expected. Both theoretical considerations and observations of the 
primordial spectra in the microwave sky can set limits on the ambiguity parameters in the model 
and highlight how the parameter range can change upon performing a symmetry enhancement 
from pure FRW to perturbed FRW backgrounds. Many of these results have been achieved when 
taking only inverse- volume corrections into account, but control over holonomy corrections is 
now sufficient for completing the picture at the effective-constraints level (see Section 9 and 
references therein). 

The purpose of this review was to show the different facets of LQC and to describe the 
work done in various directions. These indicate that, as we move towards studying the LQC of 
models with lesser degrees of symmetry, a more holistic understanding of all these developments 
is needed. This is essential for further improving our hold on this quantum theory of gravity. 
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